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We propose the addition of a new “soft-collinear” mode to soft collinear effective theory (SCET) 
below the usual soft scale to factorize and resum logarithms of jet radii R in jet cross sections. 

We consider exclusive 2-jet cross sections in e^e“ collisions with an energy veto A on additional 
jets. The key observation is that there are actually two pairs of energy scales whose ratio is R: the 
transverse momentum QR of the energetic particles inside jets and their total energy Q, and the 
transverse momentum AR of soft particles that are cut out of the jet cones and their energy A. The 
soft-collinear mode is necessary to factorize and resum logarithms of the latter hierarchy. We show 
how this factorization occurs in the jet thrust cross section for cone and kr-type algorithms at 0{as) 
and using the thrust cone algorithm at 0{a^s). We identify the presence of hard-collinear, in-jet soft, 
global (veto) soft, and soft-collinear modes in the jet thrust cross section. We also observe here that 
the in-jet soft modes measured with thrust are actually the “csoft” modes of the theory SCET+. 

We dub the new theory with both csoft and soft-collinear modes “SCET++”. We go on to explain 
the relation between the “unmeasured” jet function appearing in total exclusive jet cross sections 
and the hard-collinear and csoft functions in measured jet thrust cross sections. We do not resum 
logs that are non-global in origin, arising from the ratio of the scales of soft radiation whose thrust 
is measured at Qt jR and of the soft-collinear radiation at 2AR. Their resummation would require 
the introduction of additional operators beyond those we consider here. The steps we outline here 
are a necessary part of summing logs of R that are global in nature and have not been factorized 
and resummed beyond leading-log level previously. 


I. INTRODUCTION 


Hadronic jets play a crucial role in testing the dynam¬ 
ics of the strong interaction described by Quantum Chro¬ 
modynamics (QCD), probing the constituents of protons 
and other hadrons, and searching for signatures of physics 
beyond the Standard Model. Reliable predictions of jet 
cross sections require advanced tools in QCD, such as 
factorization of hard, collinear, and soft dynamics [T], 
the resummation of large logarithms of separated energy 
scale ratios in perturbation theory [2] , and the accounting 
of nonperturbative effects, identifying universal proper¬ 
ties when possible [3]- 

Many jet observables have been resummed to very 
high accuracy, especially with the advent of effective 
field theory (EFT) techniques applicable to jets with soft 
collinear effective theory (SCET) [H-^. Two-jet event 
shapes in e“''e“ collisions [5] such as thrust, hemisphere 
jet masses, and C-parameter have now been resummed to 
next-to-next-to-next-to-leading logarithmic (N^LL) ac¬ 
curacy HDHS], with a plethora of other observables now 
resummed to NNLL accuracy. The resummation of jet 
observables dependent on the definition of an algorithm 
have also been studied, such as angularity jet shapes m- 
m , jet masses [IMl and other jet substructure observ¬ 
ables EH [22]. One class of logarithms that have not yet 
been fully resummed are logs of the angular size R of jets 
in the definition of these algorithms. Fixed-order depen¬ 
dence on R in resummed cross sections in another vari¬ 
able has been computed, but a generic procedure to re¬ 
sum logs of R themselves has not yet been given. Ref. [23] 
has recently applied generating functional methods to 


resum leading logs of R in QCD in several jet observ¬ 
ables. The logs of R in jet cross sections and the jet 
shape Al{r/R) EHUi] were partially resummed in a soft 
approximation in |26] . and the resummation of the jet 
shape was also explored recently using EFT techniques 

m- 

In this paper we explore steps toward an EFT-based 
method to factorize and resum logs of jet radii in jet 
thrust and exclusive jet cross sections. We consider 
specifically two-jet thrust and total cross sections in e“''e“ 
collisions. For example, with the Sterman-Weinberg (S- 
W) jet algorithm [21], one obtains the total 2-jet cross 
section to 0{as), 
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where 5 is the half-angle of each cone and /3 is the frac¬ 
tion of total final-state energy lying outside the cones. 
Because S-W cones are not fixed to be exactly back-to- 
back, at 0{as) this is equivalent to the fc^-type recom¬ 
bination algorithms [21H33]- For a fixed cone algorithm 
{e.g., EH) with cones of half-angle r and energy veto A 
on particles outside, one obtains instead 
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where R = tan(r/2). In both Eqs. ([^ and <§ , (To is 
the Born cross section. Below we will identify /3 = A/Q 
and S = R in discussing the two classes of algorithms. 
Meanwhile, for a measurement of jet thrust, defined in 
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Eq. ( pj^ , the 0{as) cross section integrated from 0 to r 
asCp 


is given by 


-a,^'«(r) = l + 

(Jq ZlT 


-81ni?ln — - 1 
Qt 


—21n^r —Slnr (3) 

2Aaf^-{T) , 


for 0 < r < for fixed cone algorithm and 0 < r < 
i?V4forS-WorkT. Acr^f^s- is a nonsingular contribution 
given in Eqs. ( [4^ and ( [43| ). 

It has been a challenge to resum the logarithms in cross 
sections like Eqs. 0,© , and ([^ because it has not been 
clear how to relate the logs of i? to a hierarchy of energy 
scales in an EFT. Numerous groups have now applied 
SCET to the study of jet cross sections dependent on a 
jet algorithm, e.g., [H [13 [25H1S]- None of these have 
yet achieved a complete factorization of the logs of R in 
Eqs. 0 and 0 , as they make use of a single soft mode to 
describe the radiation with energy at the scale A = QjS. 
Ref. [13 did introduce the notion of refactorizing the soft 
function for a cross section with multiple measurements 
on jets (such as different jet masses along with the jet 
veto) but stopped short of identifying different modes 
for soft particles carrying the jet veto energy. 

In this paper we show that the logs of R can only 
be factorized completely if an additional mode is added 
to SCET below the ordinary soft scale, namely, a soft- 
collinear mode, with soft energy but collinear in its direc¬ 
tionality, with light-cone momenta Psc = {n ■ p,n ■ p,p±) 
scaling as ^ A(l,i?^,i?) or A(i?^,I,i?). We show the 
appearance of this mode in the one- and two-loop con¬ 
tributions to the jet thrust and total 2-jet cross sections 
and how it leads to more complete separation of the scale 
dependence at these orders, thus making possible resum¬ 
mation of a larger set of logs of R than before. 

We do not deal here with the factorization and resum¬ 
mation of logs that are non-global in nature |4T], that is, 
from correlated soft emissions into regions measured with 
disparate scales, e.g., the jet thrust Qt and the veto en¬ 
ergy A. The scale R also enters these ratios, as found in 
[371IM1I11I13- Recent breakthrough progress has been 
made in factorizing and resumming these observables in 
the framework of SCET [21] [^ [11] through the realiza¬ 
tion that to factorize and resum non-global logs (NGLs), 
more differential measurements of subjets, from which 
soft radiation can be emitted, must first be made. New 
modes and scales associated with the boundary of the 
full jet appear and play a crucial role in the resumma¬ 
tion of the NGLs. Similar ideas have just recently been 
outlined in [33 in the context of some of the same two-jet 
cross sections we consider here, also identifying the soft- 
collinear mode (with the less mellifluous name “coft”) 
and the necessity of additional operators accounting for 
multiple eikonal directions within jets that can emit this 
radiation. 

Resummation of the NGLs and thus a full resumma¬ 
tion of all logs containing dependence on the algorithm 
parameters A and R, and a full realization of the power 


of the soft-collinear mode, will indeed require treatment 
of emissions from multiple collinear sources within jets. 
Such a treatment, being pursued by several groups, is 
one we leave out of the scope of the present paper. We 
deal here with all the issues of implementation of the 
soft-collinear mode at the first step of the full treatment, 
namely, with one eikonal source within the full jets. We 
identify those logs that can newly be factorized and re¬ 
summed at this level with the introduction of this mode. 
We study both the jet thrust and total 2-jet cross sec¬ 
tions. In the process we tie together several dangling 
ideas in previous literature: 

• The soft mode for radiation inside jets measured 
with a thrust-like variable is in fact the collinear- 
soft (csoft) mode (not to be confused with soft- 
collinear!) introduced in [33 in the context of 
multi-jet cross sections in which two of the jets ap¬ 
proach each other in angle, increasing the virtual¬ 
ity of the exchanged soft radiation due to increased 
collinearity. This theory of SGET plus csoft modes 
was called SCET_(.. We point out here that the 
same thing happens to soft radiation contributing 
to thrust when confined to a cone. The increase in 
this soft scale by a factor 1 /R was noted in [33 ; the 
identification of the in-jet soft mode with the csoft 
mode of [33 is new. The theory of SGET+plus our 
soft-collinear modes in order to factor logs of the 
ratio between A and AR is then called SCET++. 

• The “unmeasured” jet function introduced in m 
[33 for collinear radiation in a jet not measured 
with any other variable (as occurs in total jet cross 
sections) can be constructed from a convolution of 
the measured jet function and the in-jet measured 
soft function in the jet thrust or other measured 
jet cross section. In this paper we determine that 
hard-collinear modes at the scale Qy/r and the in¬ 
jet measured soft modes at the scale Qt/R become 
one and the same mode when r is integrated up to 
its maximum value in the jet cone, R^ for a cone al¬ 
gorithm. This explains an observation in [39j about 
how to reshuffle terms between the 0{as) measured 
jet and soft functions so that they integrate up to 
the unmeasured jet function. Our observations gen¬ 
eralize all of these previous computations of unmea¬ 
sured jet functions at 0{as) and procedures relat¬ 
ing them to measured jet functions to 0{a1) and 
beyond. This can be viewed as the “0th” step in a 
procedure to resum logs in jet cross sections by first 
identifying subjets with some kind of substructure 
measurement as in E31I11I11]- 

• Previously it has been noted that logs of R in soft 
functions dependent on a jet veto [31[33ISZII3H1I13 
appear to be proportional to the cusp anoma¬ 
lous dimension [32] • While the circumstantial ev¬ 
idence from fixed-order computations was com¬ 
pelling, justification was lacking. By the introduc¬ 
tion of the soft-collinear mode, we find that both 
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the global soft and soft-collinear functions contain 
a rcuspln(Ai/(2A)) and Tcusp ln(At/(2Ai?)) pieces 
which combine to give Fcusplni?. Because the in¬ 
dividual anomalous dimensions are In ^-dependent 
we can conclude that their coefficient, and thus that 
of the Ini? in the unfactorized veto soft function, 
really is the cusp anomalous dimension, establish¬ 
ing the validity of the earlier conjectures. 

• We identify and resum all those logs of R not as¬ 
sociated with non-global radiation, distinguishing 
those coming from the hard-collinear scale QR and 
the soft-collinear scale 2AR. This is a necessary 
step to deal with global logs of R before resumming 
also the NGLs. 


• The physical origin of the argument, Qt/{2AR'^), 
of the NGLs in the jet thrust cross section now 
becomes clear as the ratio of the csoft scale of mea¬ 
sured soft radiation Qt/R and the scale of soft- 
collinear radiation 2AR, the only modes that are 
sensitive the the jet boundary. 


The necessity of additional modes such as csoft and 
soft-collinear modes can be traced to the measurement 
of multiple observables that set the collinearity of jets 
or softness of radiation emitted therefrom, in this case, 
the jet thrust r, the jet radius R, and soft veto energy 
A. Some modes are sensitive to only one of these mea¬ 
surements, such as the hard-collinear modes to r. Others 
are sensitive to more than one measurement at the same 
time, such as csoft to r and R, or soft-collinear to A and 
R. The emergence of modes simultaneously sensitive to 
multiple measurements has been explored in studies of 
multi-differential cross sections, e.g., [nmsiiis!. Hard- 
collinear modes sensitive to R were introduced in [niiTS] . 
Modes sensitive both to a soft veto A and to R, however, 
are introduced here for the first time, as is the connection 
between csoft modes and R. 

Particularly for recombination algorithms such as the 
k-r-type algorithms, jet cross sections begin to exhibit 
“clustering logs” [lH HU [50] at even in the 

Abelian part, that spoil predictions of higher-order logs 
based on renormalization-group (RG) evolution of lower- 
order logs. Thus the generality of our proposed ideas will 
depend on the choice of algorithm. We leave a discussion 
of these sorts of logs out of the scope of our treatment 
here. At 0(0^) we will stick to the thrust cone algorithm 
as in EZlISHlllS], which exhibits good factorization prop¬ 
erties at least to 0{al). As we have emphasized, the 
introduction of the soft-collinear mode is a necessary, 
though not always sufficient, step towards resummation 
of logs of R. We will find it already illuminates several 
outstanding issues, especially those outlined above. 

In Sec. |TT| we will explore the appearance of multiple 
soft scales at 0{as) in jet cross sections and introduce the 
soft-collinear mode. In Sec. Ill we will dissect the phase 
space in the 0{as) 2-jet cross sections and construct the 
separate global soft and soft-collinear functions at this or¬ 
der. In Sec. IV we will review the 0{as) jet thrust cross 


section and integrate it to obtain the 2-jet cross section, 
leading to insight into how to relate factorization theo¬ 
rems for the two cross sections to higher orders. In this 
context the complete EFT of SGET++ makes its first 
appearance. In Sec. [V] we will construct the form that 
the 0(ag) soft function for the jet thrust cross section 
must take if it factorizes into csoft, soft and soft-collinear 
pieces as we propose, and compare this form to its ex¬ 
plicit computation in [^. We find the latter result does 
indeed factor into pieces dependent on the global soft 
and soft-collinear scales, and we extract the anomalous 
dimensions of the individual csoft, soft and sof t-collinear 
factors to 0{a‘l) and even 0{a^). In Sec. VI we obtain 
the O(a^) jet thrust cross section implied by our refactor¬ 
ized formula, and relate the measured jet and in-jet mea¬ 
sured soft functions to the unmeasured jet function, ob¬ 
taining its form and anomalous dimension to 0(af). We 
also obtain the terms in the 2-jet rate that can be com¬ 
puted from integrating the jet thrust distribution with 
one collinear source of soft radiation per jet and compare 
to the recent prediction of [] 


we com- 


In Sec. |Vn 

pare some of the terms in the 0(ag) 2-jet cross section 
with the numerical predictions of the program EERAD3 
m- In Sec. jVIlIj we present a formula for the resummed 
jet thrust distribution using the ingredients we extracted 
earlier. In Sec. m we collect in summary form the main 
new results of the paper for easy reference. In Sec. 
we conclude. In the Appendices we provide a number of 
technical formulae and known results that we need in the 
main part of our discussion. 


II. HIERARCHY OF HIERARCHIES 

The jet and jet thrust cross sections in Eqs. 0, 0, 
and 0 differ in one important respect from something 
like the global thrust or hemisphere jet mass distribu¬ 
tion, which exhibit double and single logs of the sin¬ 
gle measurement parameter r or mjQ. For thrust and 
jet mass, both collinear and soft radiation are measured 
with the same parameter, e.g., t, and r acts as a phys¬ 
ical regulator on both the soft and collinear divergences 
in QGD amplitudes, generating logs of r. Factorizing 
logs of r, then, involves identifying which ones come 
from collinear divergences and which ones from soft. In¬ 
stead, in the jet rates Eqs. 0 and 0, the collinear and 
soft divergences are regulated by different parameters, 
in principle completely independent from each other. 
The soft divergences are regulated by the parameter A, 
while collinear divergences are regulated by R. These 
parameters control the amount of phase space around 
the divergent regions contributing to the cross section. 
Thus, while collinear and soft modes relevant for a global 
thrust or jet mass cross section are connected by r, e.g., 
Pn ^ Q(I, T, ffr) and Ps Q{t,t,t), giving the usual re¬ 
lation iff = Qiis amongst hard, collinear, and soft scales, 
the collinear and soft modes for the jet rates Eqs. 0 and 
0 are in principle decoupled in their scaling. 







4 


Normally, to describe a jet cross section in an EFT, we 
integrate out of QCD hard modes of momenta 


Eq. 0 are proportional to the cusp anomalous dimen¬ 
sion in Eq. (D4) according to: 


Ph ~ {Q,Q,Q) 


(4) r// = -2rcusp, rjun = rcusp, rg = 2r°^gp. (8) 


and match at a scale = Q onto a theory (SCET) of 
energetic collinear modes and soft modes. For 2-jet cross 
sections Eqs. Q and 0, these modes have light-cone 
momenta scaling as 

Pn ~ Q{1,R‘^,R) 

SCET: pf,r^ Q{R‘^,1,R) (5) 

Ps ~ (A, A, A). 


Trying to factor a jet thrust or two-jet cross-section like 
Eqs. 0,0, and 0 into jet and soft functions using this 
separation of scales yields a soft function containing logs 
of multiple scales that cannot be minimized by a single 
choice of soft factorization scale ps [HEg. From m, 
one obtains the prediction 


Here H = \C\^ is the hard function, built from the Wil¬ 
son coefficient C appearing in the matching of the QCD 
current Jqcd = q"/^q on the SCET 2-jet operator O 2 = 
Xn7±Xn, through the relation (Jqcd) = C'2(p)(02}(p), 
with the expectation value computed in some suitable 
external state overlapping with the operator [SI EU |S3] . 
The functions containing effects of the hard collinear 
radiation inside jets were dubbed “unmeasured jet func¬ 
tions” in [HEg to distinguish them from more usual 
jet functions in which some property like the jet mass 
is probed, e.g., |54[|5g. The soft function (f^e “ 
indicates integrated, or cumulative, distribution in en¬ 
ergy Eout outside the cones) contains the effects of soft 
radiation from jets outside the cones in the veto region. 

At one loop, the individual hard jSS], jet, and soft [ig 
functions in Eq. 0 are given by, 

H(Q^, ;r) = 1 + ^ (r?, In^ ^ In ^ + c\j) 

(7a) 

{QR. m) = 1 + ^ (r!} +75 In^ +c5"5f•) 

(7b) 

^5gto(A, i?, M) = 1 + ^ (r°s lni?ln ^ + 4) , (7c) 


where as = as{p), and T^p, 7 p are the 1-loop coefficients 
of the cusp and non-cusp parts of the anomalous di¬ 
mensions, and Cp are constants. The cusp parts of the 
anomalous dimensions are proportional to the universal 
cusp anomalous dimension to all orders in as, which has 
the expansion given in Eq. (D4) with coefficients given 
up to three loops by Eq. (D5). The cusp coefficients in 


In the form Eq. ( |7c| ), we do not actually know that the 
coefficient Tg of the double log is proportional to the cusp 
anomalous dimension beyond 0{as)'. it multiplies only a 
single log of p and thus technically is part of the non-cusp 
anomalous dimension in the form Eq. (7c). It has previ¬ 
ously been conjectured to be proportional to the actual 
Tcusp to all orders Eg, and confirmed to be so to 0{a1) 
by explicit computation Isa Eg. We will show below 
that the introduction of the soft-collinear mode enables 
its identification as the cusp anomalous dimension. 

The non-cusp parts of the anomalous dimensions are 
given by the expansion in Eq. (D4), and have coefficients 
for the hard and jet functions given by Eqs. (D6) and 
(D7). Meanwhile, the constant depends on the par¬ 
ticular algorithm chosen and was computed for cone and 
fer/Sterman-Weinberg algorithms in |15j . 


+r = (7 + 61n2-5^)c,, 


The hard and soft constants are given by 


/ 77 r 2 s Ott^ 

Cp = (^— 16^ Cf, 4 = 3 “^^' 


In Eqs. (7a) and (7b), it is evident that the depen¬ 


dence on the hard scale Q and the jet scale QR have 
been properly isolated, and logs in the fixed-order hard 
or jet functions can be minimized by choosing p = Q ot 
p = QR, respectively. In the soft function in Eq. (7c), 


it may appear that the scale p = 2 AVr could be cho¬ 
sen to do the same (as suggested in Eg). However, this 
choice, as we will see, does not succeed in minimizing logs 
at 0{a1) and beyond, and the anomalous dimension of 


in Eq. (7c) still contains an undesired large Ini?: 


IS = 


rflnS'veto a 


dlnp 


= ^r^lni?., 

27r 


( 11 ) 


to 0{as). As we will show, there is actually dependence 
on two scales, 2A and 2AR, which must be separated out 
to factor the soft function into pieces each depen¬ 

dent on a single scale. 

We introduce a mode in the soft sector, with collinear 
scaling with respect to the soft energy A: 


P.--A(l,l,l) 

p:,^A{1,R^,R) (12) 

p%^A{R^,l,R) 


This mimics the scaling of the hard-collinear modes in the 
original SCET Eq. 0. Once the hard and hard-collinear 
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modes are integrated out, one is left with a theory of soft 
modes, obeying the full QCD Lagrangian mm , so it is 
as if one is starting over from the full theory, but at a 
lower energy scale and one where soft gluons are sourced 
by Wilson lines instead of physical quarks and gluons. 
Within this soft QCD, one repeats the matching onto 
collinear modes with respect to the lower energy scale. 
The second theory is a “soft SCET,” and we name the 
new modes soft-collinear modes. 

Again, the soft-collinear modes in Eq. are to be 
distinguished from the “csoft” modes of |46j . which are 
soft modes whose virtuality is higher than the global soft 
scale because the small light-cone component is fixed to 
be the same as the soft while the other momentum com¬ 
ponents are forced to have collinear scaling. In [15] this 
happened when two collinear jets whose combined invari¬ 
ant mass is measured grow close to one another in angle, 
causing the exchanged soft radiation to inherit a degree 
of collinearity. Below, we will see they also arise when the 
soft radiation from a jet whose small light-cone momen¬ 
tum is measured (with thrust or mass) is confined to be 
inside a cone. In contrast, in Eq. (12), the soft-collinear 
radiation lives at a scale below the global soft scale, be¬ 
cause its large light-cone component of momentum is 
fixed by measurement of the energy /c° = -|- A:“)/2 to 

be equal to the soft scale A, while the other components 
scale collinearly. 

The Lagrangian for the soft-collinear modes looks iden¬ 
tical to that for the hard-collinear modes in the first 
SCET. They do in principle couple to a set of softer 
soft modes, e.g., with scaling in all com¬ 

ponents, which can be decoupled from the sc modes by 
the BPS field redefinition |7| , but such modes do not con¬ 
tribute anything to the jet cross sections we consider (we 
introduce no measurement sensitive to this scale). 


III. PHASE SPACE AND EFT MODES 


In this section we consider the scaling of phase space 
constraints for the cone and fer/Sterman-Weinberg algo¬ 
rithms at 0{as) and match the full QCD calculation onto 
a theory of the collinear, soft, and soft-collinear modes 
defined above. These constraints and their collinear and 
soft scaling for several algorithms were extensively stud¬ 
ied in |35|. 

The two-jet rate in QCD can be computed from 

^ 2 -jet (A, R) = ^Y. (^I r |0) I' ©alg. {X) 

^ (13) 

X {2TTf5'^{q-px), 


For a final state with hadrons i of momentum ki, 


Qcone{X) = ^ [0{9iii < r) -I- 9{9in > TT - r) (14) 

iGX 

-I- 9{r < 9in <TT - r)9{k° < A)] , 


where 9^ is the angle between the particle i and the 
cone axis ri (determined, e.g., by the thrust axis). The 
two theta functions allows particles to have any energy 
if inside the cone of radius r, and to be limited by the 
jet veto energy A if outside. For the kx or Sterman- 
Weinberg algorithms, the theta functions are more com¬ 
plicated, having to account for different possible recom¬ 
binations of particles. We will write their specific form 
for a definite number of collinear or soft particles below. 

To 0{as), the cross section Eq. (13) receives contri¬ 
butions from purely virtual, hard graphs with a qq final 
state, and from the collinear and soft gluon limits of the 
qqg final states, up to corrections suppressed by R or 

A/Q, 


a2-jet(A, R) ^(gd) + 2j(i) + .g(i)). (15) 

(To dTT 

ij(l) 

comes from purely virtual diagrams for qq final 
states, and corresponds to the squared modulus of the 
one-loop part of the matching coefficient of onto the 
SCET two-jet operator Olf |8l|52j|53]. The result is given 
by Eq. ( |7a| ). 

and S'(i) come from the collinear and soft limit 
of the QCD amplitudes and phase space constraints for 
one-gluon emission. These amplitudes and phase space 
constraints are given, e.g., in m- They are given by 
matrix elements of collinear and soft operators in SCET. 
For example, the soft function in Eq. is given by [15] 

^U(A, A,m) = ^TtJ2\{Xs\T[Y„Y^] |O)|"0:ig,(AQ , 

^ Xs 

(16) 

where Yn,n are Wilson lines of soft gluons m, and 0®ig 
are the algorithm-imposed constraints in the soft limit for 
each final state Xg. We give its 1-particle form below. 

The phase space constraints on a qqg final state in the 
limit of a gluon collinear to the n-collinear quark is 


©cone = 9 [r^- 


Pt_ 

Pg 


X 


R^ - 


pA_ 

Pg 


(17) 


and similarly for 0"one- We recall R = tan(r/2). After 
integrating over the momentum conserving delta func¬ 
tions in Eq. (13), these theta functions translate into a 


single set of constraints on the gluon momentum. 


where = q'j^q is the quark current, is the lep- 
tonic part of the amplitude (see, e.g., [56|[57]), px is the 
momentum of the final state and q is the e“''e“ total mo¬ 
mentum, in the CM frame, q = {Q, 0), and 0aig. imposes 
the phase space constraints for the chosen jet algorithm. 


= (18) 

^ Pg ^ ^ [Q-Pg?^ 

and similarly for 0”one' For the Sterman-Weinberg or 
kr algorithms, the collinear gluon phase space is given 
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FIG. 1: Phase space for cone and kT/Sterman-Weinberg jets. The phase space for a collinear (blue) or soft (red) gluon emitted 
from a quark and antiquark in the computation of the 0{as) 2-jet rate in the cone and kT-type/Sterman-Weinberg algorithms 
are shown. 


Pg jPt=Pjn 



FIG. 2: Soft phase space for one gluon. The phase space for one soft gluon emission is the same for both cone and kr- 
type/Sterman-Weinberg algorithms. The original soft phase space on the left covers the region outside both jets where radiation 
is vetoed by the energy cut Eg < A. This region is actually sensitive to two distinct physical scales. On the right-hand side, 
this region is re-expressed in terms of region sensitive to one physical scale at a time. The purely soft (or “global”) region 
covers all angles and is sensitive to the scale 2A. The “soft-collinear” regions cover gluons of energy A within the jet cones of 
angle R and are sensitive to the scale 2AR. 


by [la, 


P,n _ nf p2 Q Pg 


)■ 


(19) 


and similarly for These regions are illustrated in 

blue in Fig.and give the contribution in Eq. (151, 
given by Eq. ( |7b| ). 

Technically, the collinear gluon phase space also in¬ 
cludes the region outside the jet cones/regions, where its 
energy would be capped by the veto energy A. However, 
this double counts the region of phase space covered by 
the soft gluon, and must be subtracted out [SB] . The con¬ 
tribution of this region of the collinear phase space was 
shown to cancel against this subtraction in m- Thus it is 
not drawn in Fig.[^ In addition, such a soft/zero-bin sub¬ 
traction must also be made from the blue regions; how¬ 
ever, these give scaleless integrals in dimensional regular¬ 
ization (DR), and so we do not explicitly include them. 
They are necessary, however, in properly interpreting any 
1 /e divergences in DR as being of IR or UV origin (in ad¬ 
dition to virtual diagram contributions) [3S1 |3U [SSI [SH] . 


For jet functions with finite R and an additional mea¬ 
surement such as jet mass, the zero-bin subtractions are 
nonzero even in DR and are essential to obtain correct 
results [isjisn]. We refer the reader to these references 
for appropriate discussion. 

In the limit that the gluon is soft, ^ A Q, 
the emitting quark or antiquark determines the jet axis 
and thus automatically lies inside the jet, and these con¬ 
straints reduce to 


0 ® 


= e{e,g > r)e{e,g > r)e[pl < A) ( 20 ) 


Technically there is also a term for soft gluons inside the 
jets with no energy constraint, but in the soft limit the 
integrals over the soft amplitude are scaleless and zero in 
DR. This region is illustrated in red in Fig. [^and gives 
the contribution in Eq. ([7c|. 


Going back to the soft phase space given by Eq. (20) 
and examining Eq. (7c) and Fig.[^ we see that this phase 
space is still sensitive to two distinct physical scales. We 
separate these out in Fig. The first scale is, of course, 
the soft jet veto energy A. The other, however, is the 
soft-collinear scale AR. Although the cross section does 
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not contain any particles of this energy scale inside the 
jet cones, it is still sensitive to such modes because we 
have taken them out of the phase space. This is evident 
in Fig. and is made transparent by rewriting the theta 
function in Eq. ( [20| ): 

0™ne = Oip° < A)[l - 0{e,g < r) - 0{e,g < r )]. (21) 

Each of the three terms, represented by the Sg and scn,n 
regions in Fig.j^is now sensitive to a single physical scale. 
The contributions sum to give 


, ( 22 ) 


corresponding to the three terms in Eq. (21). 


Here we can give explicit results for the 0{as) contri¬ 
butions in Eq. (22). We will consider them at 0{a1) in 
Sec. |V| The full veto-dependent soft function in Eq. (28) 
was computed in Ref. m from matrix elements of the 


Wilson lines in Eq. (16), integrating over the real gluon 
phase space in Eq. (21). The computation there was or¬ 


ganized in such a way that we can easily read off the 
results corresponding to the individual global soft and 
soft-collinear pieces in Eq. (22), whose all-orders opera¬ 


tor definitions we give in Eqs. (26) and (27). Working in 


d = 4 — 2e dimensions to regulate the divergences, in the 
MS scheme, we find that 


^(l),bare _ 

agCp / 

/i \2Y 1 TT^X 

(23a) 

TT r(i — e) V 

2A) 6 ) 

^(l),bare _ 

agCp eT'®" 

/ /r \2<^ 1 

(23b) 

to 

1—‘ 

1 

rt\ 

\2Ar) £2 ■ 


Arguments similar to those in Refs. [S^ |31] can be used 
to identify these poles as being of UV origin. Each com¬ 
ponent can be renormalized separately, and we find for 
the renormalized functions 


where 


SW =T°,ln^^+cl 

s i, 2A 

=r°,ln2^+c 

sc sc 2Ai? 


1 

sc 1 


^pO pO 

2^5 ^ cusp ’ 






(24a) 

(24b) 


(25) 


where are Wilson lines of the soft gluons scaling as 
Ps in Eq. (12), and Xg are final states in the soft mode 


Hilbert space. The phase space constraint Eq. (20) for 


these soft modes reduces to unity since they cannot re¬ 
solve the small angle R. The global soft function Eq. (26) 


with a measurement of the total energy is equivalent to 
the timelike soft function that appears in Drell-Yan, e.g., 

in [5T] . 

Meanwhile, the soft-collinear function is given by 


1 


i?, ^) = _ ^ I r[Y-r-] 






n-Pi 


O{0in > r)] , 




(27) 


where 


are Wilson lines of soft-collinear gluons in 
Eq. (12), Xgc are final states in the (n) soft-collinear 
Hilbert space, and similarly for the n soft-collinear func¬ 
tion (which is equal to S'"^). For these modes, the soft- 
collinear gluons in the n direction only resolve the angle 
with respect to the n jet axis; the phase space is the 
complement of the sc„ triangle in Fig. Similarly the 
h-soft-collinear modes resolve only the angle with respect 
to the n jet axis. Integrating a soft-collinear gluon over 
all angles gives a scaleless integral in DR, so we know 
the contribution of one soft-collinear gluon in the region 
0iii > r is minus the contribution over 0^ < r, hence the 


minus signs in Eq. (21). 


The total veto soft function is given by 


J dEjdEgdEr,dE^5[E-Eg-E^-Eii) 

X Sg{Eg, p)S:,iE„, R, p)SUEn, R, P) , (28) 


which in this form is cumulative (integrated) in E (up to 
A). The division of this soft function into contributions 
from global soft radiation and soft radiation confined to 
the cones, and their corresponding scale dependence, was 
anticipated in the breakup of its calculation at 0{ag) in 
m- However the identification of each piece as the soft 
function for a different mode and RG running to two 
different scales 2A and 2AR was not made there. 


IV. INTEGRATING THE ONE-LOOP JET 
THRUST DISTRIBUTION 


Now the dependence of the cross section on the separate 
scales 2A and 2AR has been properly separated. 

The phase space division in Eq. (21) suggests the all¬ 
orders generalizations of the purely soft, or global soft, 
function and soft-collinear functions for veto-scale radia¬ 
tion: 


Sg{E,p) = S(^E - ^ E,), 

(26) 


Although the one-loop cross sections Eqs. 0 and 0 
are well known, we can cross check them by integrating 
the jet thrust distribution. This exercise will also lend 
us insight into a way to cross check the two-loop predic¬ 
tion for the jet rate from our calculations against known 
results for the two-loop jet thrust distribution, which we 
will carry out in Sec. |V] 

The jet thrust distribution identifies two jets with an 
algorithm just as for the total 2-jet rates Eqs. 0 and 
([^, but makes an additional differential measurement of 



























the jet thrust, one of the jet angularity shapes defined in 
m- The jet thrust is defined by 


r(J) = 


1 


2Ej 


IIIptI 




(29) 


ieJ 


where J is a jet found by the chosen algorithm, the sum is 
over particles i within the jet, Ej is the total jet energy, 
and the transverse momentum and (pseudo)rapidity 
rji are measured with respect to the jet axis. We will 
consider both the double-differential distribution in the 
thrust of the two separate jets, ti 2 , and the single dif¬ 
ferential distribution in the total thrust t = ti +T 2 . Fol¬ 
lowing the derivation in m, we find the leading power 
contributions to the jet thrust distribution can be com¬ 
puted from the factorized cross section. 


= E[{Q‘^,fi) J dtndtndkndkn (30) 


1 da{A, R) 

ao dTidT2 

n \ jr f kf\ 

X {tn, R, R, fJ-)S{kn, kn, A, R, fi) , 


X d (ti - 


where are algorithm-dependent jet shape functions 
[13 [SD] and S is the jet shape soft function, equal to at 
least 0{as) for cone and k^-type algorithms. The total 
jet thrust distribution is then given by 


1 da{A, R) 
(To dr 


= — [ dTidr2(5(r—Ti—T2) 

CTo J 


d(7{A, R) 

dri dT2 


(31) 


The measurement of jet thrust in Eqs. (30) and (31) 


induces sensitivity to a different set of collinear and soft 
modes than those in Eqs. ([^ and (12). The hard modes 
are integrated out the same way to give the hard function 
H in Eq. (30), but the EET we match onto is SCET with 


an extra set of “collinear-soft” (csoft) modes, a theory 
that was dubbed “SCET+” [i5] - 


Pen ^ Q{T,l,Vr), Pen ~ Q(l,r, Vr), 

SCET+: 

Ps ~ A(l,l,l). 


(32) 


The csoft scale Qt/R can be identified from the explicit 
computation of the soft function, given to 0{as) in [T5] 
and below in Eq. (39). Physically, it arises because the 


measured soft radiation with small light-cone component 
~ Qt inside the jet cone is being confined to angular re¬ 
gion of size i?, increasing its collinearity and virtuality— 
in the global thrust cross section, the soft scale would 
just be Qt. Confining such radiation to a cone then re¬ 
quires the large light-cone component to be ^ Qt/R? and 
p_L ~ Qt/R, as in Eq. ( |3^ . The rescaling of the mea¬ 
sured soft scale by 1/i? was identified in [15], and here 
we assert further that they are in fact the csoft modes in 
[46| . with an overall soft energy ^ Qt but with relative 
collinear scaling of the components in Eq. (32). In [l^ 


csoft modes arise because soft radiation is exchanged be¬ 
tween two collinear jets whose angular separation grows 
small. Here they arise because the measured soft radi¬ 
ation from one jet is itself being confined to a smaller 
cone. 

The jet thrust cross section exhibits dependence, of 
course, on the (veto) soft and soft-collinear scales in 
Eq. (12). We will deal with the refactorization of the soft 


function into pieces dependent on one of these scales at a 
time in the next section. The whole hierarchy of relevant 
scales is illustrated in Fig. [^ The complete EFT with 
hard-collinear, csoft, global (veto) soft, and soft-collinear 
modes we dub SCET-|__|_: 


SCET 


+-I- 


Pen 

- Q(r, 1, v^), 

^ / T T \ 

Pen ^ * 

' 

Pcs " 


or Qy- 


-A(l,l,l), 


Psc 

A{F?, 1, R) or 

A(1,A‘ 


r > 
Rj 


) ■ (33) 

To one-loop order, the jet functions J))'®' = in 

Eq. (30) are given by 


(t„, i?, p) = {tn,p)+ A J"'®- (t„, i?), 

where is the usual inclusive jet function [^ [ 

4 \9{t)\n{t/ 


(34) 


(35) 


r 3 

\e{t)p^^ 

4 r 

i p? 

t 

+ p^ L 


t/p? 


J + 


where the plus distributions are defined in App. [^ and 
Ajai®- is an additional contribution dependent on the 
algorithm. These were computed in [13 [S3 for cone al¬ 
gorithms and m for k^/Sterman-Weinberg algorithms, 
with the result: 


AJ““(t,A) = 


OtsCp 


e{t)e{Q‘^R^ -1) 


6 


47r 

0{t-Q^R^) 

I 


t + Q^R^ 


'41n 


Q2i?2 


-f 3 


(36) 


and, 


Aj'^^{t,R) = 


47r ( 


611 -I- 41n^ 


t 

1 — Xi 


t 


Xi Q2R2/] 

t 


(37) 


41n 




+ 3 


where 


Pi = 511 - 


1 - 


4t 

Q2i?2 


(38) 


Strictly to leading order in t/{QR)‘^, that is, r <C in 
Eq. (30), the algorithm-dependent corrections Eqs. (36) 
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Hard scale 


Jet scale 


Csoft scale 


Global soft 
(veto) scale 


Soft-col I inear 
scale 



iQ,Q,Q) 



Vt) 



1 

R 



(A, A, A) 




FIG. 3: Scales in SCET++ for the jet thrust cross section. The scaling of the light-cone components of momentum 


P±) 


for each mode is shown. Hard virtual modes of scale Q are integrated out. The jet scale is the same as the global thrust 
distribution; the usual soft scale is increased by 1/R due to the restriction of measured soft radiation to a cone of radius R, 
like the csoft mode of [46]. The soft veto on the energy A of additi ona l jets induces a global soft veto mode that cannot resolve 
the angle R as well as the soft-collinear modes that can, see Eq. (12l. The csoft scale here could also be below the soft veto 
and/or soft-collinear scales. 


and (37) are power suppressed relative to 

II5 


Thus 

the jet thrust cross section is 


in that limit, as in 
independent of the algorithm to 0{as). However, to ob¬ 
tain the total 2-jet rate later, we will integrate the cross 
section up to kinematically maximum allowed value of 
Tmax = R'^ (cone) or r^ax = 7?^/4 (ky), where is 

no longer power suppressed. Thus we keep the contribu¬ 
tions of both and A in what follows. 


The soft function in Eq. (30), meanwhile, receives con¬ 


tributions from soft gluons emitted within the jets and 
contributing to the jet thrust, and those outside and 
encountering the jet veto A but not contributing to r. 
Putting together these different contributions at 0{as) 
from m for two measured jets in the final state, we ob¬ 
tain 


Sik 


h- 

n ) '^n •) 


A, i?, /r) = 5{k) 

^CtsCp 


1 + 


-E 


1 

IJlR 


Att 
0{ki)^R 
ki 


4 In i?In 


4A2i? 


In 


nR 


(39) 


-I + 


which in this form is differential in kn^n but cumulative 
(integrated up to A) in the energy veto. Note that the 
parts of AJ®'^®- in Eqs. (36) and (37) for t 


cancel the total contribution of S in Eq. ( |30| ) 

above r = Tmax [H] ■ 


Putting together the hard function Eq. (7aI, the jet 


function Eq. (34) including the contributions Eqs. (36) 


and (|37|), and the soft function Eq. (39), the prediction of 


Eqs. p0[ ) and (31) for the (total) jet thrust cross section 
to 0 ( 0 ^), presented in integrated form, is 


f^c‘® (^) = — / dr 
^0 Jo 


1 ^^,da{A,R) 

0 dr' 


(40) 


= 6l(r)0(Tmax - T)<j 1 -I- 

-81ni71n^ - 1 


OisCp 


Qt 

+ B{t - Tinax)o-e^®'(Tinax) 


27r 

2Aaf^\T) 


— 21n^ T — 31nT 


where Acr))^®' is given by 


Act))'®- (t) = 6>(rmax -t) dt A J®-^®' (t, R). 
Jo 


(41) 


As noted above, the contribution of AJ‘^'® (r) above 
r = TVnax cancels against the sum of the inclusive jet 
function and the soft function contributions, so the in¬ 
tegrated distribution plateaus at its constant value at 
r = Tmax- The final term crc'®'('^max) in Eq. (40) is then 
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simply the expression in braces for r < Tmax evaluated at 
T = Tmax- For the cone and Ft/S-W algorithms, the inte¬ 
grals of the algorithm-dependent contributions Eqs. (36) 
and (371 are given by 


CisCp T + E? 


Aa““(r) = d(i?^-r)^ 61 n 


47r 


i?2 


and 




1- A 1- 


4r 


(42) 


(43) 






At the maximum values r = Tmaxi these simplify to 
2^ OisCp , 


= R^) = ^=^6 In 2 
47r 


(44) 


(r 

OisCp / 

V 4 / 

47r V 


27r2 


6 - — -61n2-l-81n^2) . 


Thus the predictions of Eq. (40) at r = r^ax are 


a™“(r = i?2) = 1-h ^|;^f-81ni?ln^ -61ni? 

ZTT \ Q 


-H61n2- l) , 


(45) 


and 


R 




= — = 1 


asCp 

27r 


2A 


f — 8 In i? In —— 6 In i? 
\ 


5 - 


27^ 

3 


(46) 


agreeing with Eqs. 0 and 0. 

Note that the differential jet functions J‘^^^ {t,R,fi) 
when integrated up to t = Q^T^ax do not by them¬ 
selves re prod uce the unmeasured jet functions J{R,p) 
in Eq. (7b I —the coefficient of the double log dif¬ 


fers. However, after combining R, fj) with 

S{kn, kn, A, R) through Eq. (30), the total r distribution 
integrated up to r = Tmax does equal the total 2-jet rate 
Eq. 0 or Eq. 0. This is because the wide-angle radi¬ 
ation in a measured jet with jet thrust r belongs in the 
csoft sector in Eq. ( [30| — it cannot become very energetic 
and preserve a small value of r for the jet—while in an un¬ 
measured jet, hard collinear radiation can go all the way 
up to the angle R (or R/2) and remain in the allowed 
collinear jet phase space. This difference was noted in 
m- Thus some of the logs appear in a different sector— 
csoft or hard-collinear—depending on whether r is mea¬ 
sured or not, but their total contribution to the 2-jet 
cross section remains the same. The transition between 
the two descriptions is relevant to formulating a complete 
all-orders form of the factorization theorem that resums 


all logs, global and non-global, in the 2-jet cross section 

(c/. [illEH m]). 

A slightly different approach was taken in |39j . where 
a definition of the measured (or, there, “unintegrated”) 
jet function was adopted so that it does integrate up to 
the unmeasured (or, “integrated”) jet function. This was 
done by adopting the power counting i? ~ t from the be¬ 
ginning, whereas [H] began with t R. Ref. [53] effec¬ 
tively combined the measured jet function J0^s (t„, i?,/i) 
in Eq. (34) with the part of the soft function Eq. (39) de¬ 


pending on the in-jet measured radiation with momenta 
kn into a single object that by itself integrates up to . 
Their soft function then retains no dependence on the 
measurement r and is the same for cross sections with 
measured or unmeasured jets. We have shown above 
that keeping the scales associated with r and R sepa¬ 
rate {t R from the beginnning but then integrating 
the total cross section up to r ~ R) does in fact repro¬ 
duce the total 2-jet cross section to leading power in R. 
Both approaches lead to the correct fixed-order total 2- 
jet cross section. However, in the approach of [39], one 
could not study the resummation of logs of the jet thrust 
T itself, although as a path to the total 2-jet cross sec¬ 
tion it is essentially equivalent to ours. The /i-dependent 
logs in the jet and soft functions of [39] are only those 
of fi/{QR) or /i/(2A) (no soft-collinear scale was identi¬ 
fied), with no /^-dependent logs containing the jet mass 
p2 or soft momentum k relevant to the resummation of 
logs of T remaining. Eor the goal of simply obtaining the 
total fixed-order two-jet rate, however, their approach 
succeeds in organizing the pieces of the jet thrust cross 
section in a way that can be integrated simply to give the 
total two-jet rate. We will explore the generalization of 
the relation between measured and u nmeas ured jet/soft 
functions to 0(a^ 


and higher in Sec. VIB 


V. TWO-LOOP JET THRUST SOFT FUNCTION 

In this section we explore the factorization of the two- 
loop soft function for the jet thrust cross section in 
Eq. (30) into csoft, global (veto) soft, and soft-collinear 


functions. The total jet thrust soft function was calcu¬ 
lated to 0{al) in [38]. Ideally, one should compute the 
individual factors from their operator definitions, e.g., 
Eqs. ([2^ and ([^, the results of which to 0(as) we 


gave in Eq. (23), and verify that they reproduce the re¬ 
sult of [38]. We do not do the EFT computation explic¬ 
itly at 0{a^) in this paper, leaving it for future work. 
Rather, we show that the result in [35] does indeed fac¬ 
tor into pieces sensitive individually to csoft, global soft, 
and soft-collinear scales, and that these agree with the 
generic form these functions should take if they satisfy 
the appropriate RG evolution equations (RGEs). There 
are leftover non-global pieces sensitive to multiple scales, 
and we show these remaining pieces agree with known 
coefficients for leading and subleading NGLs at 0(ag). 
From 1381 we can then extract the anomalous dimensions 
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of the individual global factors to 0{a'j.). 

Beyond this, we give generic arguments for how the 
csoft, global soft, and soft-collinear anomalous dimen¬ 
sions should be related to each other and to other soft 
functions in the literature (such as those for hemisphere 
thrust and Drell-Yan) to all orders in as, and show 
our extracted anomalous dimensions indeed, nontrivially, 
satisfy these relations. The generic relations allow us to 
extract the anomalous dimensions of all the soft functions 
even to 0{a^). These results are already highly nontriv¬ 
ial and are a strong consistency check of our refactoriza¬ 
tion framework. From an explicit EFT computation of 
the csoft, global soft, and soft-collinear functions from 
their operator definitions such as Eqs. (26) and (27), one 


should be able to confirm our extracted anomalous di¬ 
mensions from their UV divergent behavior and repro¬ 
duce the results for the renormalized functions we give 
in this section. 


A. Refactorization and RG Evolution 


In the soft function for jet thrust in Eq. (39), there is 


dependence both on the csoft scale of measured in-jet soft 
radiation k/Ras well as the scales 2A and 2AR associated 
with the energy veto outside the jets, illustrated in Ei g.[3| 
We notice the dependence on the jet veto scale in Eq. (39) 


arguments of the ratios of the various soft scales [HI 03] . 
They come from correlated emissions into separated re¬ 
gions of phase space. itself is a convolution of two 
identical pieces.) For the soft function for the total jet 
thrust cross section in Eq. (31), 


S(^k, A, R, |x) — / dkjidkn^i^k kji kji^S{kji, kji, A, R, |x). 

(51) 


We will distinguish the two soft functions Eqs. (50) and 


(51) by the number of their arguments. The refactoriza- 

from the jet veto-dependent 
and the separation 


tion of the A:-dependent Si 
parts was already proposed in [15 


of Sng and identification of the coefficient of the leading 
(double) NGLs was discussed in 03]. The identification 
of the subleading (single) NGLs in this soft function was 
made in [38]. The separation of the soft and soft-collinear 
functions at the scales 2A and 2AR here is new. The full 
set of scales we have in the jet thrust cross section is 
illustrated in Eig. [^ 

We will not provide here a formal derivation of 
Eqs. (50) and ([^. However, because an explicit two- 


loop calculation of the total jet thrust soft function in 
Eq. (51) is available |38| . we can check if it does in fact 
obey the form Eq. ( |51[ ) in a consistent manner up to 
0{a'^), and proceed to compute the prediction arising 
from Eq. (50) and the full factorized form of the cross 


2 -jet cross section, and thus separates into soft and soft 
collinear contributions as in Eqs. (22) and (24). In fact. 


is the same as in the soft function Eq. {Jc^] for the total section Eq. (30) for the resummation of (global) logs to 


the whole 1-loop jet thrust soft function Eq. (39) can be 
separated into three categories of terms: 

Sikn, kn,A, R, /r) = 6{k) [1 + 5(1) (A, ^l) + 2Sil'> m)] 

+ (47) 


where were given in Eq. (24), and where 5Y’' is 
given by m, 


Ai) 


all orders in Us at NLL and NNLL accuracy. 

It will be most straightforward to work with the 
Laplace transform of Eq. (51) in both k and E. We form 
the double transform: 




(52) 


5(i/, A, i?,/r) = / dke-''^ dEe' 

Jo Jo 

X Si^{k/R,^if®Ss{E,iJL)®Ssc{ER,^f®Sng{k,E,R ), 

which disentangles the convolutions, leaving us with 
S{v, A, R, y) = Sf^ivR, m)'5s(A, h)S^^{X/R, ii)Sng , (53) 


S[^{k/R,y)='^{5{k)c^^-^-^ 


Ak 


2Fn 


e{k)\a{k/{iJLR)) 


^R [ k/(fJ-R) 


(48) 


where 


cL = -^Cj.. 


(49) 


The form Eq. (47) suggests a generalization to an all¬ 
orders form of a refactorized jet thrust soft function, 


Si^kji , kn , A, 7?, /x) — / dE Sini^kn/ ll) iSin {J^nj Ri fJ) 

Jo 

X Ss{E, fj,) (g) Ssc{ER, fl) (g) Sng{kn,n, E, R) , (50) 


where the (g are convolutions in the jet veto energy mea¬ 
surement variables and/or in kn^n- Sng contains the left¬ 
over fixed-order NGLs, which begin at 0{a^), and are 


where 


Sin{ER,tl) = 

poo 

/ dke-''’^Sin{k/R,ti) 

Jo 


Ss{X,h) = 

pOO 

/ dEe-^^Ss{E,^j) 

Jo 

(54) 

Ssc{X/R,n) = 

poo 

/ dEe-^^Ssc{ER,fi) 

Jo 



These functions should satisfy ordinary RGEs of a mul¬ 
tiplicative form, 

fl) = Jin{fJ^)Sin{l^R, fJ.) 

= 7ss(m)>^s(A,/x) (55) 

fjj— Ssci^XR^ fjj^ — 7sc(/^)*5'sc(A-R, /i), 
dji 
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where the anomalous dimensions take the forms 
7in(Ai) = -2rcusp[Q:s(Ai)] In —-h 7 in[as(^)] 

(qj 

= 4rcusp[as(Ai)] In -h 7^*[os( a*)] 

lsc{p) = -2rcusp[as(Ai)] In + 'ysc[asi^J -)], 

where the 7 [as] terms on the right-hand sides are the non¬ 
cusp parts of the anomalous dimensions, and the propor¬ 
tionalities of the cusp parts to Fcusp are deduced from 


the one-loop results Eqs. (251 and (48). 


Among the objects in the jet thrust cross section 
Eq. (|30|) , the anomalous dimensions of the hard function 
mm and the jet function (which is the same as 
for the inclusive jet function) [65] as well as the cusp 
anomalous dimension m issj are all known to 0{al). 
(We give these in App. |^) The jet thrust soft function 
for the thrust cone algorithm, and thus its anomalous 
dimension, were computed to 0(al) in |37l |38|. Thus 
far, whether it factorizes as Eq. (50) and what the forms 
or the anomalous dimensions of each individual piece 
are remain unknown. Under the assumption that they 


obey Eqs. (55) and (56), we can extract these individual 


anomalous dimensions from the form of the soft function 


B. Two-Loop Soft Function 


A function that obeys the multiplicative RGE Eq. (55) 


with anomalous dimensions of the form Eq. ( [56| ) whose 
expansions in Us take the form Eq. (D4) must take the 


form given by Eqs. (Dl) and (D2) (see, e.g., l67]l. 

that is. 


a. 


Si{xi,fj,) = 1 + — 


\n^{gLXie' 




(s)‘ 




+ (ri+S^rO) ln2(/ix,e^®) 

+ (7*^ + ‘^c\j3o) ln(/rx,eT'®) -f 

where we have already used that the 0{as) non-cusp 
anomalous dimensions 7 ° = 0 for each function in 


Eq. (54). The arguments Xi for each function are: 


Xiri h'R , 


A 

X q G - i 

ss 2 




2R 


(58) 


The cusp parts of the anomalous dimensions in Eq. (57) 
are determined by 


T- — r — —r 

-L in — -L sn — A (7 


r — 2 r 

■L S.9 - C 


(59) 


and the one-loop Laplace-space constants (see Eq. 
are given by 


cj = cj 


T'^ — 


(60) 


where c] for i = Sin was given in Eq. (49), and for i = S. 


and Sec in Eq. (25). 

Since the cusp anomalous dimension appearing in 
Eqs. (|5^ and (Dl) is known to two (in fact, three) loops, 


and the non-cusp anomalous dimensions and constants 
to one loop, up to NLL(') resummation (cf. [TTJ [S7]) of 
(global) logs is already possible. To go to NNLL accu¬ 
racy, the non-cusp anomalous dimensions are needed to 
two loops. For NNLL' accuracy and higher, the two-loop 
constants are also needed (besides calculation and resum¬ 
mation of NGLs |H|). In this section, we will extract the 
two-loop non-cusp anomalous dimensions 7 !, from 
known results for the full (but not fully refactorized) jet 
thrust soft function in the literature [35] • 


Using the expansions Eqs. (Dl) and (D2) for Si 


fixed-order expansion of S{v,\^ R, gi) in Eq. (53) up to 
0{a1) can be written 


in,ss,sc 


( 57 ) and multiplying them together according to Eq. (50), the 


5'(:/, A, i?, At) = 1 -f — 2ro f- In^ URve'^’^) + In' 
47r L V 


AtAe'’'^ 


-In' 


2 AtAe'’'® TT 


2R 


- y 1 + 2c(„ 




1„2 AT! _ 1„^ 


/xAe' 


7 e \ 2 


+ 2ro/3o(-ln^(At-Rt^e'''®) + In' 


2 
3 


2cl 


2R 

In^ A, i?At) 

2 AtAe'>'® 


2 2R 

2 . 


+ 2 Ti -f To ^2cin -I- -I- 2c^^ ~ y ) (~ ) + In 


— In' 


AtAe' 


7E 


AtAe'^^ 


2R 

AtAe'>'® 


+ 2(7in + 2/3oCin) ln{pRve^‘^) + ( 7 ]^ -I- 2/3oc^s) In - h 2 ( 7 ]^ -|- 2/3oc],n) In 


(61) 


■ Co 


Our goal is now to extract the (thus far) unknown two-loop anomalous dimension coefficients 71 ^ We will not 
extract the two-loop constant c| in this paper. 
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To compare Eq. (61) to the expression for 5'(/c, A, i?,/x) given in [38) . we perform the inverse Laplace transform of 
Eq. (61) back to momentum space, using the formulae in Eq. (B5), 


S,{k,k,R,^l)=L~^CT4 


1 
v\ 


S{iy,X,R,fi) 


(62) 


for the cumulative soft function in k and A, resulting in: 

,2 . ^2 

3 

,3 M 


S^k,A,R,fi) = l + ^^ 




o/3o ( 


-In' 


3 k-R , ,3 k 


In —— In' 

2A 2KR 

fj,R 


-SilHk,A,R,fi) 


(63) 


^R 


+ 2 (ri - ToyCf) (- In^ ^ + Ini? In^) - ^(ro)2(ln2 ^ + In^ i?) 

k 


/ii? 


+ 2(7in + 2/3oCin - 8C3ro) In -^ + (tL + 2/3ocL) In ^ + 2 ( 7 ]^ + 2/3ocl^) In , 

where ^ is the constant term in this momentum-space total soft function. 

We can separate the purely Abelian (agCp and a'^Cp) terms and non-Abelian {CfCa, CfTfUf) terms of Eq. (63), 


S^{k, A, R, ^Ji) = Sc, {k, A, R, k) + {^) ^S^niik, A, i?, ^i ), 


(64) 


where 


Scp{k,A,R,fj,) = 1 + — . 

Att L 


'(■ 


2 fJ-R 

k 


2To ( - In^ ^ -I- In i? In 




iA^Rj 

2 


iAV 


4A2i?y 


+ (S) {2(ro)"(-ln"^+lni?ln 
+ 2ro^c^(in^ X - ^" 4 ^) - X + X + ’ 


(65) 


and 


S^^lik, A, i?, /i) = ^ro/3o (- In^ X + X ■ 21/?) + X + ^ 


^R 


+ 2(7in + 2/3oc(„) In ^ + ( 7 L + 2^ocL) In + 2 ( 7 ]^ -h 2/3ociJ In . 


(66) 


2A 


2Ai? 


In [38], an expression is given for the non-Abelian terms of the 2-loop soft function both for arbitrary R and also in 
the limit that R —>■ 0. We are interested only in these terms that do not vanish as i? —>■ 0. These terms can be directly 
compared to Eq. ( | 6 ^ for extraction of the unknown anomalous dimensions 7 !^ We quote the formula from [35] 
App. [^ In terms of known cusp anomalous dimension and beta functions coefficients, their result can be reorganized 
into the form: 


® (- f + - '”X) + 4®;) 


-(2) 

^raA[38] 


-h 


-h 


-h 


( 56 C 3 - 
r/1616 


^R 

X 

1616 

~w 


LV 27 
( 56 C 3 


448 47r2 

CfCa + c^CfTfuj + —PoCf 

- -^CfTfuj - ^PoCf 


In 


pR 


56(3^CfCa 


1616 

XT 


448 


47r^ 


CfCa + -^CFTFUf -\— —PqCf 




16. 

y 


I6C3 — 7) CfCa + —CFTFUf -\ — —PqCf 


2A 

In yy+ 

^ 2AR^ 

SX 

3 


In 


(67) 


2Ai?2 ■ 


On the last line we have separated out those logs which sions into two separated phase space regions. These cor- 
are non-global in origin, coming from correlated emis- 
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( 2 ) 

respond to the Sng term in Eq. ( 661 : 

2 ^ 




16, 


2AR‘^ 

8tt^ 


( 68 ) 


Cf (^-16C3 “ 3 )*^"^ + ~ 3 ~^° 


In 


2 Ai?2 ’ 


the leading coefficient of which was computed in [43] , and 
the single log coefficient in [371|3S] . (Ref. |33|, however, 
only identified one power of R in the argument of the 
NGLs). The computation of Eq. (671 i n 1371138] together 


with the factorization conjecture Eq. (50) confirm that 
the argument of the NGL is, in fact, k/{2AR^) [43], which 
we notice is the ratio of the measured in-jet soft scale 
k/R and the soft-collinear scale 2AR in Fig.[^ From the 
“in-out” and “in-in” NGL coefficients computed in [45] 
and the results of (33], we can also form the correspond¬ 
ing non-global contribution to the double-differential jet 
thrust soft function in Eq. (50), 




Cf[{-1GC3-^)Ca 


16 


Tpn/ 


2AR2 
8^2 


X (In 


2AR2 


+ ln 


2AR2y 


V 


(69) 


In the limit i? —)■ 0, no NGL of kn/kn appears [43], as the 
phase space for soft gluons inside the two cones vanishes. 
A resummation of the NGLs requires a more advanced 
factorization theorem using technology such as that in 

[mill]. 


Comparing Eqs. (66) and (67) and using Eqs. (25) and 


(|49[) for the one-loop constants we are able to 


read off the non-cusp anomalous dimensions: 


iL = -2"fin = -27ic 


= Cf 


/1616 
V 27 


56C3)c.-^T.n,-^ 


(70) 

The relation amongst the three anomalous dimensions in 
Eq. (70) should, in fact, generalize to all orders in as- In 


the i? = I limit, the in-jet regions each become a whole 
hemisphere, so Eq. (56) directly implies for the non-cusp 


anomalous dimensions: 


Tin — Themi 
Tss “ ‘^^sc 


(71) 


We can derive several more useful relations amongst 
anomalous dimensions. To start, jss is the same anoma¬ 
lous dimension as for the timelike soft function that arises 
in threshold resummation in Drell-Yan, a property that 
was in fact already noted in This follows from the 
definition of Ss{E,fi) in Eq. (26). The only difference 
with Drell-Yan is the direction of the Wilson lines (in¬ 
coming vs. outgoing), which does not affect the anoma¬ 
lous dimension (nor the value of the function itself to at 


least O(a^) [55]) and so 

+ ‘^Iqq + 7 ss = 0 , (72) 

where jqq is the non-cusp part of DGLAP evolution of 
quark parton distribution functions, which for example 
appears in the Altarelli-Parisi splitting function as 

Pqq{z) = + y,Ja,]5(l - z) + ... , (73) 

I-*- - z)+ 

where the ellipses denote terms that are non-singular as 
z —>■ 1. Then, we use the condition for the consistency of 
factorization in thrust (e. 3 ., [SUET] El HO]), 

7 i/ -I- 27 j -I- 27hemi = 0 (74) 

and the relation satisfied by "fH,j,qq in DIS as a: —>■ 1 [65j . 

7H+lJ+lqq = 0, (75) 

which together imply 


TJ ^yhemi — T^g ■ 


Finally, Eqs. (72), (74), and (76) together imply 


^Themi Tss — 0 ■ 


(77) 


Combining this with Eq. 0 , we derive the all-orders 
relations among non-cusp anomalous dimensions, 


7hemi — 7in 


= 7sc = - 


"fss 

2 


(78) 


We give 7hemi to three loops using the known results 
quoted in Eq. ( |D8[ ). 

The relations Eqs. 0 and ( |74| imply the consistency 
of anomalous dimensions in our jet thrust factorization 
theorem Eqs. (30) and ([^, 


JH + 27j -|- 27in -|- 7ss -|- 27 sc — 0 , 


(79) 


which in this form is satisfied by both the cusp and non¬ 
cusp parts. 


VI. TWO-LOOP FIXED-ORDER CROSS 
SECTIONS 


With the results for the two-loop soft functions in 
Sec. [V] and the known results for the two-loop hard and 
jet functions in App. [^ we can construct our prediction 
for the logarithmic terms in the full jet thrust cross sec¬ 
tion to two loops, and integrate it up to r = Tmax to 
obtain our prediction for the logarithmic terms of the 
two-loop jet rate for the thrust cone algorithm. 
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A. Jet Thrust Cross Section 


^ng ij^n i kni E ^ FC) . 


(80) 


We have the ingredients to build the jet thrust cross 
section Eq. (30) differential in each jet’s thrust ri ^2 or 
the total thrust r in Eq. We will begin by working 

with the double distribution Eq. (30), though it is easier 


to express it in terms of its integral over ri_ 2 - Using the 
refactorization of the soft function Eq. ([50|) , we have that 


(Jc{ti,T2\E,R)= f dT[ f 

Jo Jo CTo aTidT 2 

X jdtndkn9(Tl-^-^'^J^^^■{tn,R,^i)Sin(^^,^J^ 
X jdtndknO ("^2 - ^ (tn , R, M)'S'in ^ 

X J dEsdEsc9{A -Eg - Esc)Sb{Es, fi)Ssc{EscR, 


We will use this factorization to resum logs in the pertur¬ 
bative expansion of the cross section in Sec. |VIIl1 Eor the 
purpose of computing the fixed-order cross section up to 
0{a^), it is convenient to package pieces of the factoriza¬ 
tion formula Eq. (80) together, those contributing to ti _2 
and those depending on the soft veto A. 

Eirst we compute the convolution of and S’in- We 
will compute the product of their Laplace transforms and 
then inverse transform back to momentum space. The 
Laplace tr ansf orms of J{t, R, /i) and Sin{k, R, /i) take the 
form Eq. (Dl) to 0{a^). The jet function also includes 
an algorithm-dependent component AJ’^*®’. 

Up to 0{a'^), the Laplace t rans form of th e measured 
jet function is given by Eqs. (Dl) and (D2) with Tj = 
^Tcusp and jj = 2 : 


poo 

j :'®- {u,R,^m)= dt j :'®- {t, i?, fi) 

Jo 


/o 

= 1 + ^ 
Att 


— (^In^ ^ 




6)2 Q2 


+ cj + (^2, R) 




4 fJ. I'e ' 


-k 

-k 


Ti 

2 


7g + 

2 2 


Q2 ^ 




In' 


3 




(" 


2/3o / 1 To TT 
- [cj + —— + AJ, 


2 6 

( 1 ) 


In 


2 




2 6 




In^+a^. + AjS, 


(81) 




y,R) 


All of the coefficients here are known (see App. except for the two-loop algorithm-dependent correction . It 

is not expected to contribute to the singular part of the r distribution (or later to logs of R in the integrated 2-jet 
cross section) at O^a"^), just as it does not at one loop. The cross terms between the one-loop AJ,^jg and logs from 
the one-loop jet and soft functions, however, do. So we keep them. We will not compute the constant Cj in this paper, 
nor compute explicitly the Laplace transform Abut keep it in abstract form until we return to momentum space. 

Meanwhile the Laplace transform of the measured soft (csoft) function Sin is given to 0{ag) by 


Sin{iyR,fi)= / dke-'''^/^Sin{k/R,fi) 

Jo 

2 iiRiye^^ 


= 1-k 


47r - 


-Tn In" 


Q 

ri + To 


TT 

-Tn 


+ Ci, 


6 


In 


-(s)‘ 

2 fiRve'^^ 


(ro)2 4 ^iRve!'^ 2 3 


(82) 


2/3o(c-n 


-Tn 


6 


In 


nRve'^’^ 

Q 


+ Cs, 


where all the coefficients are also known except the constant cf^, which we do not compute in this paper. The two-loop 
non-cusp anomalous dimension 7 ^,^ is now known through our extraction in Eq. (70) from the results of [38] . 

Multiplying the Laplace transforms Eqs. (81) and (82), we use Eqs. (B5), (B 6 ), and (B7| to inverse transform the 


logs and products of logs of ly. The result of this exercise is 




= 1+^ 

47r 


Tnl^ln^ 


•( 
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Yilr,2 

Q 


QtJ 
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-k Cj -k cfn + (t, R) 




QtJ g 2 
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+ [ri + ro(cj + cl„) I ( — In^ -In' 


1,2 1 '" 

2 "" QV *“ QtJ 


1 , „1 M 1 ^ 1„2 M_, 2 tJ-R \ 


(75) 


0\2 


■IjPo 


1 , 2 

4 Q'‘t 


°V12 QV 6 Qt 6 QR 
2 7j + c5(7j + 2/3o) + c?n7j- 


In 


7in + ^clj/So — 2C3ro + ro7j^ 


+ ^ dr'AjW (r-r',i?)|ro 


In 


Qr 


^2 


cl - Tn 


12 


) + const. + A (t, i?) 


(83) 



rin(QV7^2)i 

Q 

9ro ^ 

'\ii{Qt'/{^ lR)y 

75 + Wo Q^ 

r 

.)) 


Cp-T' 

+ Vr 

. QT'/ifj-R) 

+ 2 /i2 

[qW'\ 


The algorithm-dependent one-loop contribution AJ^|g^ is given by the integrals of Eqs. (36) and (37). For the cone 
algorithm, 


AJ,^W(Ai?) = 6C'fln(l + ^). 


The convolutions on the last line of Eq. (83) are given by Eq. (BIO), with the result 

\n{Q^T' / ^^)'\ Q 


dT'AJ^ll{T-T',R)l^r, 




= scp.; 7 






- 2 rn 


J + 

,2 


' iiR 


Ri{Qt' |{^lR)) 


Q^t 
2 


-Tr 


2 Qt 


In^ Inf 1 -I- 

fiR V R-^J 


Qt 


Qt'/{^ iR) 
2 Lis 


-I + 


2 ^JR [qV' 


„ , — 2 In Li 2 f--I- 2 LIt . „ , 

i?2y ^R ^\t + r:^) ^\r + R^)\ 


(t^j?) 

(t^) 


7±^ri„«7i„(i + i_)_Li. 

2 L ^JR \ R^J 




(84) 


(85) 


Meanwhile, the A-dependent parts Ss 0 in Eq. (80) can be combined into the function, 

2 




' 1=1 + 
,2 4 


2 ro In i? In -h 2 c^^ 


+ (S) + 3!^°^° ^ 2177 + 

2 

+ 2[Ti+ ro(ci, -h 2clJ] lnR\n^^ + Wo{c\s + 2c^J ^ -f ( 7 ^ “ 4/3ocL) Ini? -f consts.j , 


where cj 


were given in Eq. (25), and we have used obey an RGE with the non-cusp anomalous dimension 


Eq. (70) to relate 7 , 


g,, to each other. We note that 5'^eto 
almost takes the form of product of two multiplicatively 
renormalized functions given by Eqs. (D1 1 and (D2), ex¬ 
cept for the ( 47 r 2 / 3 )(ro)^ In^ i? term, which arises be¬ 
cause S^eto is actually built out of a convolution of pieces 
in the energy variable E, integrated up to A. 

We note that 5'veto has no non-cusp anomalous dimen¬ 
sion even at two loops, since from Eq. (70), 


7veto — 7ss T 27 sc li ■ 


(87) 


Thus, without refactorizing S'veto into soft and soft- 
collinear pieces as in Eq. (50), one would not know 


that the 7^3 In R term in Eq. ( 86 ) is really built from 
the anomalous dimensions of two functions Ss sc which 


yjg = —2^1^ and thus have no way to resum the asso¬ 


ciated logs of R. Also in the form Eq. ( 86 ), one would 


not know that the coefficients of the double logs of R are 
given by the cusp anomalous dimension; this can only be 
deduced from the refactorized form Eq. (50) and the RG 


evolution of the single-scale dependent pieces Ss^sc- 
Now, combining the pieces Eqs. ( [M| an d ( 86 ) together 
with the hard function given by Eq. (|D1|) and Sng given 


by Eq. (69), we can organize the cross section Eq. (80) 


into two /r-independent factors. Expanding the cross sec¬ 
tion in as{Q), we obtain 

tTc(Ti,T 2 , A,i?) = crc(ri; A, i?)crc(T 2 ; A, i?), ( 88 ) 

where to 0 ( 0 ^) for the cone algorithm. 
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a““(r;A,i?) = 1 + 


asiQ) 

4:TT 
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ro(21iii?ln^ - 
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+ 2ro/3o(: 


2A 2 
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In^ 




- ^lnT + 4„t + 6C'Fln 




(89) 


- -In^ 
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ro75 
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2A 2AR 
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(ro)^^(in^ 
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ri + To + QCf In^l + 
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2 In i? In — -In^ —- 

2A 2 R^J 

47r^ 


+ 81n^i? 
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/ 4\ 8 47r" 

(-8C3 - ^)CfCa + -CpTFUf + —PoCf 


l) 


^ ss 

2 ' 2 
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75 
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In T 


In 


Qt 

2KR? 


&CfR(^ 

ro75g-C.rg-6&r„Li,(^) 

(c 5 + eCi. ln(l + ^ In r + 2c\^ ^ + ^ 3^0 ^ + ^cL In ^ 




+ const. 


where 


-'H 


Ctot — “ 7 “ + Cj + Cin + 


- R + 9 R ('QD'l 

2 ' '-j ' '-^in ' 2 ’ Weto *-'ss ~ •^'^sc ’ 

and where we have dropped purely non-singular terms in 
Eq. (89). We also used the consistency relation Eq. (79) 
to replace in the single log terms. 


B. Measured and Unmeasured Jet Functions 

We observed in Sec. [IV] that integrating the one-loop 
jet thrust distribution up to t = r^ax reproduces the 
total one-loop two-jet rate for the cone and kr/Sterman- 
Weinberg algorithms. In fact we can identify the “un¬ 
measured jet function” in Eq. Q with the convolution 
of the measured jet function in Eq. (30) and the part 
of the refactorized soft function in Eq. (50). Namely, we 
find that 


p^max r / 

{QR ~ J dr JdtdkS\^T — 


X J^^s-it,R,RSin{k/R,fj.), (91) 

Explicitly at 0{as), 

dr' Jdtdkd(T'--^-^^J^^s {t,R,n)Sin{k/R,fi) 


= 1-f 


where 


CXg 

47r - 




Q'^t 
2 


QtJ 


(92) 


+ Cj + Cj^jj -I- A J°ig (r, R) 


AJf 


rQ^r 

,{t,R)= / dtAR^^-{t,R). 

Jo 


(93) 


At r = Tmax, we obtain for the cone (r„iax = dR) and 
kT(Tinax = i?^/4) algorithms, 


7“ 


I 5in = 1 


Q:, 

47r - 


Eoln^ ^ ^ 


+ 7 / In 

QR QR 

+ (7 - — + 61n2jC^^ 


(94) 


and 


0 = 1 ■ 


a 
4 ft - 


T 7 1 2 M , 01 M 

1 0 In-h 7 7 In- 

“ QR QR 


(--") 


Cf 


(95) 


which equal the one-loop unmeasured jet functions given 
in Eq. (7b), with constants in Eq. (§. 

Let us continue the evaluation of the right-hand side 
of Eq. (91) to 0{al). Although an independent direct 
computation of Jun to 0{a^) does not yet exist, we can 
construct the right-hand side of Eq. (91) and check that 
it is consistent as a definition of Jun at least to 0{al). 
Namely, Jun should obey the multiplicative RGE 


7^ ^^JuniQdi^ if) — 7Jun(/^) Jun (1^77; M) J (96) 

with an anomalous dimension given by 

7Jun(/^) = 2rcusp[Q^s(/4)] In 4" 7Jun [CTs (m)] ■ (97) 


That is, Jun must take the form Eq. (Dl), which solves 
Eq. (96), with the cusp part of the anomalous dimension 
proportional to Ecusp according to Eq. ( |9^ . 

We have evaluated the right-hand side of Eq. (91) to 
0{al) for arbitrary r in Eq. (83). Setting r = Tmax = 
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there, including the contributions from A Jaig. in Eq. (85 1 , 
gives Ju™® constructed according to Eq. (91) to 0(c^. 


The relations that we need are given in Eq. (B 8 ). Using 
these, we obtain for the 2 -loop unmeasured jet function. 


Tcone _ -I j_ 


^0 7 ^ ~ ~6—6 In 2) Cf 


QR 

\2 


QR 


fas\‘^({To)^ 4 ^ / 0 , 2,o\, 3 M 


QR 


QR 


7 } + Tin + (Tj + 2/3o) (^7- - —h 6 In 2 j Cp 




L +ro ( 7 - ^ +6 In 2 ) Cf + ^ +7j/3o 


In 


(98) 

2 

QR 


In —— -I- const s. 

(cj It 


We note this is precisely the form of a function obey¬ 
ing a multiplicative RGE, as we expect from the naive 
factorization in Eq. §), which is given by the expansion 
Eq. (D1), with 


r 


Jun 


= r, 


cusp 5 


Tjuii 7j 4” Tin 


(99) 


That Jun constructed according to Eq. would sim¬ 
plify to this form obeying a simple multiplicative RGE 
was by no means obvious from the beginning, and is 
a strong consistency check of the notion of the unmea¬ 
sured jet function. The non-cusp anomalous dimension in 
Eq. (1^ can be constructed out of Eq. ([70| and Eq. (|D7| , 


H 

7Jun = 'IJ Tin = ^ ; 


( 100 ) 


where we used Eqs. (78) and ( [7^ to relate yin = —7ss/2 
This can be computed to 0{a^) 
This demonstrates 


and = - 27 j -|- 7 , 
using the results quoted in App. 
that the unmeasured jet function introduced in m does 
not simply have the same anomalous dimension as the in¬ 
clusive jet function (as was already seen from the differing 
cusp parts of the anomalous dimensions). Physically this 
is because the unmeasured jet function includes energetic 
radiation at the edges of the jet that is treated as csoft 
when measured with a thrust or mass measurement, but 
becomes of the same scaling as the hard collinear radia¬ 
tion when integrated up to r = i?^, as we see in Eq. (321. 

Once contributions from Sm are lumped together with 
J))'®- when T = to make the Jun, the veto-dependent 
parts of the soft function Ss 0 in Eq. (80) do not con¬ 
tribute anything to the non-cusp anomalous dimension 
according to Eq. (87), resulting in the simple propor¬ 
tionality between yjun and 'Jh in Eq. (100). 


In the jet thrust cross section, we can still distinguish 
between the global logs coming from RG running between 
the scales in Fig. [3l a nd the NGL in the fixed-order total 
soft function Eq. (|50[). In the integrated total 2-jet rate, 
however, this NGL, in Eq. ( 68 ), becomes a log of 2A/Q: 


In^ 


Qt t= 


R‘ 


2 AR 2 


In^-^ 

2A 


( 101 ) 


which is indistinguishable from other global logs in the 
2-jet rate. It will require a more complete dissection of 
the dynamics in the jets to fully resum both types of logs, 
such as those recently advanced in [HlllllllllS]. Our 
treatment here provides key steps that must be part of 
such a complete analysis, allowing a resummation of all 
the global logs appearing in the jet thrust cross section, a 
clearer understanding of the relation of the integrated jet 
thrust and total 2 -jet cross sections, and a determination 
of the evolution of the unmeasured jet function to two 
loops. 


C. 2-Jet Rate 

The total 2-jet rate can be computed from the double 
integral of the double differential jet thrust distribution 
Eq. (ISOl): 


rR 


<72°j“(A,R)= / dri / Jt2 


1 da{A, R) 
(To dTidT2 


( 102 ) 


for the cone algorithm. The 2-jet rate is determined by 
this double integral instead of a single integral of the total 
jet thrust distribution Eq. (311 since the size of both jets 


is limited to radius R, independent of the other jet. (The 
single integral would allow one jet to be fatter than R and 
one narrower.) From Eq. (30) and the construction of the 


unmeasured jet function in Eq. (91), we find the total 2- 


jet cross section is built from products of the unmeasured 
jet functions: 


(7®°"®(A, R) = H{Q\ fi)Jun{QR, Pif 

X ^^(A, fi) (g) SsciAR, I 


ccone 
' '^ng 


(103) 

(A,R). 


We computed the combination Sveto = Sg ® in 
Eq. ( 86 ). The non-global part is given at 0{a‘l) by 
Eq. ( 68 ) with k QR^. The total jet rate can be ob¬ 
tained from Eqs. (| 88 | and (89) with t = R?. We obtain 


^cone 1 
^2-iet = 1 


Ois{Q) 

-( 

47r 


.( 1 ) 




r( 2 ) 


(104) 
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where cr 


( 1 ) 

cone 


is 


where = (12 In 2 — 2)Cf- As for the two-loop terms, 


= 4ro In i? In ^ - 27 ^ In i? -h 


(105) 


o 9 o 2A n 9 2A / 9 2A 2A 9 \ 

cr^one = SFq In^ i? In^ -I- 8 ro 7 j In^ i? In 4ro/3o / 


Q 


Q 


- 4[ri + roC,Ve] In A In ^ + 2 {fjf + In^ R 




(l6C3 + l)CA-fTFnf-^/3o 


Q 


1 2A 

2A 


(106) 


+ [1h + - 27jCeone “ 4(cju„ J/3o] In A - 2{cl, + 2c\^)j3o In -^ + consts., 


where we used Eqs. (70) and (100) to relate the non-cusp anomalous dimensions in front of the single Ini?. Plugging 


in the explicit values of the coefficients in Eq. (|106p, we obtain 

2h 


= dCl I (321n2 + 481n + 18 - ^ ] In^ i? + 


Q 


2A 9 

(8 — 48 In 2) In -I -1- 27r^ — 24C3 — 36 In 2 

Q 2 


In R 


+ 4CfC^| (y In ^ + 11) In^ R + 


■44. 

-T 


In^ 


2A 


'47r^ 


27r" 2 2A (2 _ ^ 

-In^-h - + 4C3 - 

3 Q V3 


llTT^N, 2A 

- In — 

9 / Q 


Q 


Q 


268\, 2A 57 „ 

-) In-h I 2 C 3 - 221n2 

9 / Q 2 


Q 


Ini? 


(107) 


\ /4 

47r^ N 

2A j 

Ini? - ( - - 

- In — ^ 

/ V3 

9 / 

Q J 


consts. 


This agrees with the result given in [1^ for all the In" i? 
terms, as well as the CfTfu/ ln(2A/(3) terms, all of 
which can be accurately predicted by our factorization 
theorem Eq. (103) with soft/soft-collinear emissions from 
one collinear direction. Emissions from two or more sep¬ 
arate Wilson lines in the fundamental and adjoint rep¬ 
resentations will affect the logs of 2A/Q, beginning with 
single logs in the Cp and CfCa color factors at 0(ag), as 
found in the results of |45| . The nf term still comes from 
emissions from one Wilson line and is thus accurately 
predicted. 


VII. COMPARISON WITH FIXED-ORDER QCD 
CALCULATIONS 


We compare terms in our two-loop computation of 
the jet rate in the thrust cone algorithm to the full 
QCD fixed-order calculation implemented numerically in 
EERAD3 I^. We expect all the In i?-enhanced terms 
as well as the CfCa ln^(2A/(5) term, which is the well 
known leading NGL, and the ln(2A/(5) terms, 

which can be computed from soft emissions from a sin¬ 
gle Wilson line per jet, to be correct. To demon¬ 


strate, we subtract all the In i?-enhanced terms and 
the CfCa ln^(2A/(5) term from the numerical data and 
check whether the difference has at most single logs of 
2A(Q. We can only test those terms with an explicit 
ln(2A/(5) since we need to make a measurement of en¬ 
ergy or other variable on the final state to be able to 
easily plot the EERAD3 data. From the remainder, we 
can extract the coefficients of the ln(2A/(5) terms, whose 
accurate prediction requires additional operators in our 
factorization theorem. 

Focusing on the terms that are A-dependent, we com- 

f’2'1 

pare the differential distributions da cone/dx plotted as 
functions oi x = — InA/Q for various choices of i? 
(Fig. 1^. After subtraction the differences in da^ne/dx 
should be constant in x, which allows us to extract the 
\a(2A/Q) coefficients. We also check the i?-dependence 
by plotting with fixed A/Q (Fig. [^. To compare di¬ 
rectly with the EERAD3 output, the distributions are de¬ 
composed into three color structures proportional to N^, 
1/Nc and n/. 


da 


( 2 ) 


dx 


O' + o I /TVc On f, 


(108) 


where and Oi/j^^ are linear combinations of the Cp 
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4 5 6 7 8 9 10 


-In A/Q 


FIG. 4: Comparison between Eq. | |106| l (solid lines) and the EERAD3 output (dots) for ctn^, and a„^, plotted (upper) 

as functions of — In A/Q with fixed R. The linear dependence on — In A/Q is shown as a check for the In i?-enhanced terms 
predicted by our factorization theorem. Note the numerical sensitivity for small values of R and A/Q. The differences (dots, with 
error bars) between Eq. (1061 and the EERAD3 output are plotted (lower) as functions of — In A/Q with R = 0.32. Performing 
fits within the fit regions (red), we extract the coefficients Cvc = —76 ± 8, Ci/m^ = —1.8 ± 0.3 and Cnf ~ 10.0 ± 0.1. 





FIG. 5: Comparison between Eq. (1061 (solid lines) and the output of EERAD3 (dots) for ctn^, cri/N^ and cr„j., plotted as 
functions of — In iZ with fixed A/Q. i'he quadratic dependence on — In iZ is shown again as a check for the In iZ-enhanced terms 
predicted by our factorization theorem. 


and CfCa terms derived from Eq. (106). 
dependent terms can not be seen in dacthe 
comparison is left for future work. 


The pure R- 
/dx, and their 


Power corrections are expected to become more signif¬ 
icant when QR^ > A. Thus for smaller values of R, the 
agreement of the singular terms between QCD and SCET 
at two loops can only be seen at smaller values of A/Q. 
We evaluate the jet rate using EERAD3 for iZ = 2" x 10“^ 
with — 2 < n < 5. The calculations require the technical 
cutoff in EERAD3 to be set to 10“^®, and we have included 
six billion partonic events in the numerical integration. 
In the singular regions of R and A/Q —>• 0, the upper 
plots in Fig. |^show the linear dependence of daiohe/dx 
on — In A/Q, and the ones in Fig. show the quadratic 
dependence on — In R. These show the agreement of the 
In iZ-enhanced terms predicted by our factorization the¬ 
orem. 


The differences between our analytic predictions and 


the EERAD3 outputs are shown in the lower plots in Fig.|^ 
The differences are close to constant in a; = — ln(A/Q), 
We perform fits within the fit regions (red) using con¬ 
stant functions to extract the coefficients Cat,., and 

Cnj: of the ln(2A/Q) terms in the jet rate. We use the 
EERAD3 output for R = 0.32 in the fits, but the results do 
not change much when including more data with smaller 
iZ’s. The bands are fit uncertainties corresponding to 
the values of per degrees of freedom deviated from 
the minimum by at most 1. We obtain, 

Cn^ = -76 ± 8 
Ci/jv, =-1.8±0.3 

= 10.0 ±0.1. (109) 

The coefficients extracted from the result given in [1^ 
are 

V3 18 


21n2-121n2 2±C3) = -85.1 
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Ci,N^ = ^(l-21n2 + 61n22-5C3) =-1.56 
/'47r^ 4\ 

--j = 10.2, (110) 

where Cnj is also predicted by our refactorization the¬ 
orem. Note the good agreement of Cnj with the value 
extracted from EERAD3. The values of Cn^ and Ci/^^ 
extracted from [45] are also consistent with the ones ex¬ 
tracted from EERAD3 within the numerical uncertainties. 


VIII. RESUMMED JET THRUST CROSS 
SECTION 

The solutions of the RG equations obeyed by the hard, 
jet, and soft functions in the jet thrust cross section in 
Eq. (301, with the soft function refactorized as Eq. (501, 


give rise to the prediction for the resummed (integrated) 
jet thrust cross section (see, e.g., lezl), 

tTc(T,A,i?) = — [ 

o-Q Jo dr' 

V Q / V Qt ) 

12aJ [2Ar) H{Q ,HH)d{T)e{A) 


X j(dQ - 
X 


1 


1 Mss 
In —.^1 Mss 

2 A ’ ^ 


)"5in(5n 


2 g7Ef2 


1 Min 

r(l-O) 


Ms 


2 e.'iET 


G5„g(QT/(2Ai?2)), 

where 


2AR'^"'') r(l-T) 

( 111 ) 


Al(MffiMJ , Min 1 l^ss ; l^sc^ fi) = + 2Kj{^j,fi) 

+ 27^in(Min;M) ^ ss (Mss ; fi) -f 2Kscip-sci m) 

( 112 a) 

11 = ^{f^j7 Min; m) — ^^j(mJ; m) T 2cajj,(^j7,, /i) (112b) 

^ f^so M) ^ss(Mss; M) ~b ‘^^scif^so M) ; (112c) 


and the individual evolution kernels are given by 
= —jpKFKri^^F, fJ-) + 

wf(mf, m) = -KFWr(MF, m) ; 

where 


(113) 


K^{^iF,^J.)= [ ^rcusp[as(M')]ln-^ , (114) 

Jfe M Mf 

wr (mf , m) = / ^ rcusp [as (m' )]; 

K■yEi^J-F,^^) = [ ^7fK(m')]- 

Explicit expansions of these functions in ag can be found 
in many places, e.g., m- The values of the coefficients 


Jft^f in Eq. (112 1 for the various functions in Eq. (111) 
are given by 


Jh — I'in — Jss — jsc — 1 7 jj — 2 , 

kh — 4, sij = 2 , Fiu = 2 , sigg = 4, Aligc = 2 . 


(115) 


In Eq. (Ill I we include the NGLs coming from S'ng in 


Eq. ( | 6 ^ simply at fixe d ord er, here to 0{al). 

We have written Eq. (111) using the formalism found in 


[S3IZI] that turns the Laplace-transformed functions into 
differential operators with respect to fl, T to facilitate 
the transformation back to momentum space. With the 
results of Sec. especially Eq. (70), we have enough 
information to resum logs of r, R, and A/Q to NNLL 
accuracy (modulo the NGLs). 

In Fig. we il lustrate the integrated jet thrust cross 
section Eq. (Ill) in the thrust cone algorithm at Q = 
100 GeV, R = 0.2, and A = 10 GeV, with global log 
resummation at NLL and NNLL accuracy. At N^LL ac¬ 
curacy, the cusp anomalous dimension and running of Og 
are included to 0{ag~^^), the non-cusp anomalous dimen¬ 
sions to ©(ttg), and the fixed-order hard, jet and soft 
matrix elements to 0(ag“^). For N^LL' accuracy, the 
matrix elements would be included to 0{a^) [TTII^, but 
we do not do this in this paper. We evaluate Eq. I® 
at the central values for the scales, 

M = M// = <3 , MJ = QVr , 

Mss — 2 A, Msc — 2A.R. 

To produce the result of the unrefactorized cross section, 


e.g., with the soft function Eq. (7c) at 0{as), we set 

Mss = Msc = 2A'/R , (117) 

the choice one would make in an attempt to minimize 
logs in the 0{as) soft function (which we saw in Eq. (63) 
does not work at higher order). For this exercise, we do 
not include the NGLs Eq. ( | 6 ^ at unprimed NLL and 
NNLL accuracy, counting it as part of the fixed-order 
soft function, but of course in a full analysis we would 
have to include them. 

For the estimation of perturbative uncertainty, we per¬ 
form several sets of scale variations. First, we vary all 


scales in Eq. (116) together up and down by factors of 
2. Next, we vary each of the r-dependent scales /ryin up 
and down using the formulae: 


Mj(F;ej) = 
Min(F, 6in) — 


1 + ej 1 - 


1 + Sin ( 1 — 


(i-^)*] 0 ^ 

'■ ' nia v ' -■ 


max 

T 


^ I mav X 


. R 


(118) 


and vary eyin each between (—0.5, 0.5) {cf. [TTl 17^1. 
Recall Tmax = for this cone algorithm. Finally, we 
vary the veto-dependent scales gigs,sc each up and down 
by 2 (for the unrefactorized case, we vary the single soft 
veto scale in Eq. (117)). We then add all these variations 
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With s-c refactorization No s-c refactorization 



FIG. 6: Resummed integrated jet thrust cross section. The integrated jet thrust cross section at Q = 100 GeV, R = 0.3 
and A = 10 GeV using the thrust cone algorithm computed from Eq. (Ill I is shown with (left) and without (right) the 
refactorization of soft- and soft-collinear modes in the soft function Eq. (50l. The central values of scales on the left are 
^,as = 2A and /Csc = 2KR, and, on the right, /Css = Msc = 2Ay/R. The uncertainty bands come from the scale variations 
described in the main text. In both plots there is good convergence from NLL to NNLL accuracy (in global logs). Without 
refactorization, however, the overall scale variation is significantly larger, indicating better control of the perturbative series on 
the left. (The NGLs are not included in these plots.) 


in quadrature. These give the uncertainty bands plotted 
in Fig.§ 

For a more robust uncertainty estimation, prohle func¬ 
tions should be used for the jet and csoft scales [II1IZ3] 
instead of the canonical scales in Eq. (116 1 , but this level 


of detail is not what we are after here. We simply illus¬ 
trate the broad qualitative effect on the quality of the 
resummation of logs of r and R with and without refac¬ 
torization using the soft-collinear mode. In both plots 
in Fig. 1^ we see good convergence as uncertainty is re¬ 
duced from NLL to NNLL. In the right-hand plots, with¬ 
out refactorization, the overall uncertainties for the same 
amount of scale variation are a bit larger. This indicates 
better control over the perturbative series in the left-hand 
plot with refactorization and resummed logs of R. 

We do not include a numerical check against full QCD 
results for the jet thrust cross section here; the validity of 
the prediction of refactorization for the fixed-order global 
logs up to 0{a1) and the leading NGLs in cone and Vt- 
type algorithms was tested against EVENT2 [TIES] in 
[43j . where excellent agreement was found. 

Tmax we COuld 


By integrating Eq. (Ill) up to r = t, 


obtain the jet rate with logs of global origin resummed 
to NNLL, although in this case accounting for the NGLs 
is more essential as they also turn into logs of 2 A /(3 as 
in Eq. (101), indistinguishable from the global logs. We 


leave the proper factorization and resummation of the 
NGLs, applying, e.g., the formalism of [211 El mi , to 
future work. 


IX. SUMMARY OF NEW RESULTS 

Here we collect, for convenience, the key new results of 
our paper, up to 0{a^), in terms of anomalous dimension 


and beta function coefficients given in App. [^ All of 
these are in fact known to 0{al). 

The collinear-soft function appearing in Eq. ( [SOj ) for 
measured soft radiation in a two-jet event confined to be 
inside a cone of radius R, in integrated form from 0 to k, 
is given up to 0{ag) by 

= 1 + g(-r„k2 T+A) (119) 
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2 ^° 


/ii? 

T 


r^ln^^-Tro/Soln^ 


/ii? 

IT 


^ ' fiR 


+ (-r,-cJ,r„-)L(r„)i)i„i'^ 

+ + 2Cin/lo - ICalro)^) In G + 4n 


where and is given in Eq. (70). 


The Laplace transform of Eq. (119) is given by the form 
Eq. 


The global soft function, defined by Eq. (26), for soft 


energy outside two back-to-back jets integrated up to the 
veto A < A, is given up to 0{al) by 


.5,^(A,M) = l + ^(2roln2^+cL) 


( 120 ) 




+ ^Fo/Soln^ ^ 
° 2A 3 2 A 


477-2 

2ri + 2cLro-—(Fo)' 
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In 
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2A 


+ (tL + 2css/3o — 16C3(ro)^) 


where and 7 ]^ is given in Eq. (|70|. We leave 


undetermined. The Laplace transform of Eq. (120) is 


given by the form Eq. (57). 
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The soft-collinear function, defined by Eq. (27), for 
soft energy outside two back-to-back cones of radius R 
integrated up to the veto E < A is given up to O(a^) by 

m) = 1 + g (-ro In" ^ + cl) ( 121 ) 


+ (S) 


vT^ln^ 


2AR 2AR 


+ (-ri-cLro-y(ro)2) in^ 

+ ( 7 L + 2cL/3o - 4C3(ro)") In 


2AR 


2AR 




= {7r^/6)Cp and '-yl is given in Eq. (70). 
We leave cL undetermined. The Laplace transform of 


Eq. (121) is given by the form Eq. (57). 


The non-cusp parts of the anomalous dimensions of 
the csoft, global soft, and soft-collinear functions in 


Eqs. (119), (120), and (121) are all related to each other 


and the known anomalous dimensions of the hemisphere 
[iniEG] and time-like Drell-Yan [BT] soft functions. They 
are related to all orders by 


7hemi — 7in — 


Is 


= Is 


( 122 ) 


given explicitly to 0(al by 

-0 =0 


7s 


(123) 


71c = C'f[(^-56C3)c^ 


448 27r2 

—Tpuf - —/3o 


and to 0(al in Eq. (D8). 

The convolution of Sg 0 Sl to give the total unrefac- 
torized veto soft function, integrated up to Eout = A, is 


given to 0(al by Eq. (86). It has no non-cusp anoma¬ 


lous dimension, but has leftover large logs of R. 

The unmeasured jet function that appears in the 
total two-jet rate in Eq. (m and is constructed according 
to Eq. (foTl), is given to 0(as) for cone-type algorithms in 
Eq. (|94[) and for S-W or k^-type algorithms in Eq. (|9^ , 
and to 0(ag) for the thrust cone algorithm by Eq. (98): 


Jsr(QH, rt = I + S (r"+7 In ^ + c;.„) 

+ (r, + c^jieTo + -( 7 °)^ + 7j^o) In^ 


+ ( 7 ,kn + clnehj + ‘^M) ^ 


,2 

cone 


(124) 


where 


= (7 -I- 6 In 2 — 5Tr‘^/6)CF, and we leave 
undetermined. Jun has a non-cusp anomalous dimension 
given by 


7Jun — 7J 4” 7in 


7g 
2 ’ 


(125) 


which is given to O(a^) by Eq. (D6). 

The 0{al) fixed-order jet thrust cross section for the 
thrust cone algorithm predicted by these results is given 
in Eq. (^, and the formula that resums globa l log s 
of r and R in this cross section is given by Eq. (|111 ). 
Eq. (106) gives the terms in the fixed-order 0{al) to- 
tal two-jet rate for the thrust cone algorithm predicted 
by our formulas. All the log” i?-enhanced terms and the 
ln(2A/Q) term for the Uf color structure are predicted 
accurately by this result. The pure ln(2A/Q) terms for 
the other color structures require inclusion of additional 
subjet operators in our factorization theorem Eq. (103) 
to be completely predicted. 


X. CONCLUSIONS 

In this paper we have worked out several ramifications 
of introducing a new mode, the soft-collinear mode, into 
SCET, allowing the separation of scales 2A and 2AR cre¬ 
ated by the measurement of a soft veto energy outside 
cones of radius i? in a jet cross section, as well as the 
resummation of the logs of R due to the ratio of these 
scales. In the jet thrust cross sections, we found the pres¬ 
ence of hard-collinear, collinear-soft (csoft), global (veto) 
soft, and soft-collinear scales, and leading to the theory 
SCET_|_+ to factorize and resum logarithms of ratios of 
these scales. 

Studying the jet thrust cross section using the thrust 
cone algorithm, we extracted for the first time the two- 
and three-loop non-cusp anomalous dimensions of the 
refactorized soft functions for in-jet measured (csoft) ra¬ 
diation, global soft radiation at the veto scale, and soft- 
collinear radiation at the veto scale in the cones. 

We related the jet thrust cross section and the total 
two-jet rate, deriving a relation between the measured jet 
and in-jet soft functions and the “unmeasured jet func¬ 
tion” in the total two-jet cross section, and as a bonus 
extracted its two-loop and three-loop non-cusp anoma¬ 
lous dimensions for the first time. 

We compared the predictions of the refactorized jet 
cross section with the output of EERAD3 to 0{al. Fi¬ 
nally, we provided a formula the resummed jet thrust 
cross section (modulo NGLs) and performed global log 
resummation to NNLL accuracy. 

We should emphasize once again, as we began, that 
our treatment here was not aimed at a resummation of 
the NGLs that necessarily appear in these exclusive jet 
thrust or total n-jet cross sections. Technology to do 
this in the framework of SGET is rapidly developing, es¬ 
pecially with the recent work of mHHIllllS]. These 
developments have made clear the necessity to identify 
subjets or multiple collinear Wilson lines emitting soft 
radiation into other regions in order to resum NGLs. 

What we have explored here is the necessity of intro¬ 
duction of a soft-collinear mode to factorize and resum 
even the global logs of R arising from the ratio of soft- and 
soft-collinear scales at 2A and 2AR and various ramifica- 
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tions and insights following from this step. It is the first 
step towards a full resummation that would include the 
NGLs. Not only have we factorized and resummed the 
global logs in the jet thrust and total 2-jet cross sections, 
we tied together some loose ends and made new connec¬ 
tions among ideas in the existing literature. Namely, as 
mentioned, we clarified the relation between measured 
jet and in-jet soft functions for jet thrust with the un¬ 
measured jet function of [Til US] appearing in total jet 
cross sections. We identified the in-jet measured soft ra¬ 
diation as actually the csoft modes of [15] and its merg¬ 
ing with the hard-collinear mode in total jet cross sec¬ 
tions as r —7^ Tmax ~ R^- We established on firmer foot¬ 
ing earlier conjectures about the proportionality to the 
cusp anomalous dimension of In R terms in soft anoma¬ 
lous dimensions for jet cross sections, and extracted from 
the computations of [ST] [38| the two-loop anomalous di¬ 
mensions of the individual in-jet measured soft function, 
global veto soft function, and soft-collinear functions. We 
believe the insights here even about the global logs in jet 
cross sections provide firmer ground on which to build 
frameworks to resum all logs, global and non-global, in 
jet cross sections. 


where 


F{x) = dx'fix '). (A2) 

Thus, integrals of regular functions against plus distribu¬ 
tions obey, for a: > 0, 


dx'[9{x')f{x')]+g{x) 


= ( dx'f{x')[g{x') - g{Q)\ F g{Q)F{x), 
^0 


(A3) 


so, e.g., 


dx' 


9{F) In" x' 


X' 

r,.' 


9{x') 


r In" x' ln"+^ X 

= dx' — [gix')-gm + -^m- 

Jq ^ Tl 1 


(A4) 


Appendix B: Laplace Transforms 
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Appendix A: Plus Distributions 


The plus distribution prescription can be used to make 
a singular function into an integrable distribution. For a 
function /, its plus distribution is defined by {e.g., [73]), 


{fix)]+ = lim ^ [9{x - e)F(a;)] 

= lim - e)/(2;) - S{x - e)F{x)] , 


(Al) 


Here we collect results for the Laplace transforms and 
inverse Laplace transforms between the logs 

L = In — , L = ln(j^e^®) (Bl) 

T 

The Laplace transform is defined by 


F{v) = ^{F}{i^) = / dre-^"F(r). (B2) 

Jo 

The transforms of the logs we need in this paper are given 


by 

.if|l|=^ (B3a) 

^ ^ ly 

if|L}=^L (B3b) 

^ ^ jy 

= + (B3c) 

= 1|z 3^’^^^2C3| (B3d) 


The inverse Laplace transforms are given by 

_ p+ioo 1 _ 

.^-^{F}{t)= (B4) 

J 7 — 200 
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where 7 lies to the right of all the poles of F in the 
complex plane. The inverse transforms we need are 


- 2-k^ In^ 2 + 2 In'" 2 + 


2 IC 3 
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To derive the last three relations, we used the results for 
the convolutions in momentum space, 


.Vi 60 

Some other particular transforms we need in the main 


text that can be derived from Eq. (B5| are 
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Appendix C: Jet Thrust Soft Function 

The non-Abelian terms of the two-loop jet thrust soft 
function in Eq. (64) that survive in the limit i? —>■ 0 , 
given in [38], are 
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(We have translated the nota tion [HH] from —>• k, 

r —>■ w —?► A.) In Sec. |VB we reorganize this formula 

into a refactorized form that matches Eq. (50). 


Appendix D: Anomalous Dimensions 

The hard function and Laplace transform of the jet 
and soft functions in Eq. (30) have the expansions 




n—0 


where 


(Dl) 


so jj = 2 and jp = js = 1- The denominators Qj^s hr 
the logs depend on the precise definition of the measure¬ 
ment in the jet or soft function. 


The expansions in Eq. (D2| contain the coefficients in 


the expansions of cusp and non-cusp anomalous dimen¬ 
sions in as, given by 


OO 


Ecusplas] - r„, 7 [as] - 51(47^) 


n—0 


n—0 


The first few coefficients for the cusp anomalous dimen¬ 
sion are given by 


Eo = l 

po 

Fi = ^Ll + ^Lp+c^p 

Jp JF 


1 

+ IT 


Jp 

1 


(D2) 




Fp 


H —:—(jp + Cpl% + ‘^Cpl3o)Lp + Cp . 
jF 


Eor F = H,J, S, the logs Lp are of the form 


L// = In — , Lj = In ■ 


Q^j 


Ls = In 


IJLve " 

Qs 


(D3) 


To = ACp , 


(D5) 


El = 


4Cf [( 


r2 — 4:Cp 


r/67 
iT “ 
z245 

ivlT 

418 
55 


24 20 

Ca —^TpUf 


1347r2 

27 

407r^ 


9 

IItt^ 

' 

56(3 


22(3 \ 


3 


+ (■ 

+ {-f + 16 C 3 )CpTpnf-^Tjn} 


The non-cusp anomalous dimension coefficients we need 
in this paper are given up to three-loop order for the hard 
and jet functions by 
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The anomalous dimensions of the csoft, global soft, and soft-collinear functions are all related by Eq. (78), 7hemi = 
7 in = —7 ss/ 2 = 7sc, and can be determined through Eq. (74) or Eq. (79) to be equal to 7sc = —7 j— or, explicitly, 
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The running coupling as(/i) itself is given to three-loop 
order by 
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